NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


THESIS 

LOW  SIDELOBE  REFLECTOR  ANTENNA  DESIGN 
FOR  COMMUNICATIONS 

by 

Nelson  A.  Armas 
September,  1993 

Thesis  Advisor:  _ ^ David  C.  Jenn 

Approved  for  public  release;  distribution  is  unlimited 


Best 

Available 

Copy 


form  Approved 
OMBNo.  0704^188 


REPORT  DOCUMENTATION  PAGE 


Public  reporting  burden  tor  theis  collection  ot  intormetion  i«  estimated  to  average  1  hour  per  response  including  the  time  tor  reviewing  instructions  searching  ensting  data  sources 
gathering  and  maintaining  the  data  needed  and  completing  and  reviewing  the  collection  ot  intoimat.on  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  ot  this 
collection  ot  mtormation  including  suggetions  tor  reducing  this  burden  to  Washington  Headquarters  Sendees  Directorate  tor  tntormation  Operations  and  Reports  I2i5  Jefferson  Davis 
Highway  Suite  1204  Arlington  VA  22202-4302  end  to  the  Office  ot  Management  and  Budget  Paperwork  Reduction  Project  (0704-0 IBS'!  Washington  DC  20503 


1.  AQENCY  USE  ONLY  (Lmv«  Stenltl 

2.  REPORT  DATE 

3.  REPORT  TYPE  AND  DATES  COVERED 

September  23.  1993 

Master's  Thesis 

4.  TTTLEANDSUBTTTLE 

6.  FUNDING  NUMBERS 

LOW  SIDELOBE  REFLECTOR  ANTENNA  DESIGN  FOR 
COMMUNICATIONS 

6.  AUTHOR(S) 

Armas,  Nelson 

7.  PERFORMING  ORGANIZATION  NAMES  (S)  AND  ADDRESS(ES) 

a.  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School 

Monterey.  CA  93943-5000 

8.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADORESS(ES) 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 


Ilic  view  s  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  position  of  the  Department  of  Defense 
or  the  US  (iovemment. 


12a.  OISTRIBUTION/AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  relea.se;  distribution  is  unlimited. 

13.  ABSTRACT  (  Maximum  200  words ) 


The  currcnl  anienna  used  in  the  Marine  Radio  AN/MRC-142  is  a  focal  point  fed  reflector  with  a  first  sidclobc  of  -10 
dB  A  new  antenna  design  is  required  w  ith  a  sidelobe  reduelion  of  another  10  dB.  and  a  minimum  impact  on  the  system 
upgrade  cost  An  improved  feed  consisting  of  a  fed  dipole  and  disk  in  the  vicinity  of  a  parasitic  ring  was  designed  for 
use  witli  llie  current  reflector.  The  disk  acts  as  a  ground  plane  for  the  dipole,  with  the  ring  narrowing  the  beamwidth. 
resulting  in  a  stronger  edge  taper  on  the  reflector 

An  existing  computer  code  based  on  the  the  method  of  moments  w  as  used  to  compute  the  gain  and  the  radiation  pattern 
Circular  parabolic  reflectors  of  36  inches  and  48  inches  were  analyzed  The  results  show  that  the  36  inch  reflector 
cannot  satisfy  the  conditions  over  the  entire  frequency  band  (1350  to  1850  MHz)  simultaneously.  The  48  inch  reflector 
came  \  er>  close  to  satisfying  all  the  conditions,  although  the  sidelobe  level  is  off  by  approximately  2  dB  at  wide  angles 
in  the  H-plane 


14.  SUBJECT  TERMS 

AN/MRC-142  Marine  radio.  Low  sidelobes.  Reflectors.  Method  of  Moments 


16.  PRICE  CODE 


15.  NUMBER  OF  PAGES 

76 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 


18.  SECURITY  CLASSIFICATION 
OFTHISPAQE 


18.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


20.  LIMITATION  OF  ABSTRACT 


UL 


UNCLASSIHED 
NSN  7540-01-280-5500 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev.  2-89) 

PrescnbwJ  by  ANSI  Std  239-18  298-102 


Approved  for  public  release,  distribution  is  unlimited 


LOW  SIDELOBE  REFLECTOR  ANTENNA  DESIGN  FOR  COMMUNICATIONS 

by 

Nelson  Armas 

Lieutenant,  Ecuadorian  Navy 


Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGfNEERING 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September,  1993 


Approved  By: 


C  _ 

Dr,  David  C.^enn,  Thesis  Advisor 


W«v 


Dr  Rama  Janaswamy, 


Second  Reader 


Dr  Michael  A.  Morgan,  (^airman 
Electrical  and  Computer  Engineering  Department 


II 


ABSTRACT 


The  current  antenna  used  in  the  Marine  Radio  AN/MRC-142  is  a  focal  point  fed 
reflector  with  a  first  sidelobe  of  -1 0  dB.  A  new  antenna  design  is  required  with  a  sidelobe 
reduction  of  another  10  dB,  and  a  minimum  impact  on  the  system  upgrade  cost  An 
improved  feed  consisting  of  a  fed  dipole  and  disk  in  the  vicinity  of  a  parasitic  ring  was 
designed  for  use  with  the  current  reflector.  The  disk  acts  as  a  ground  plane  for  the  dipole, 
with  the  ring  narrowing  the  beamwidth,  resulting  in  a  stronger  edge  taper  on  the  reflector. 
An  existing  computer  code  based  on  the  the  method  of  moments  was  used  to 
compute  the  gain  and  the  radiation  pattern.  Circular  parabolic  reflectors  of  36  inches  and 
48  inches  were  analyzed.  The  results  show  that  the  36  inch  reflector  cannot  satisfy  the 
conditions  over  the  entire  frequency  band  (1350  to  1850  MHz)  simultaneously  The  48 
inch  reflector  came  very  close  to  satisfying  all  the  conditions,  although  the  sidelobe  level 
is  off  by  approximately  2  dB  at  wide  angles  in  the  H-plane 
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I.  INTRODUCTION 


A.  BACKGROUND 

An  antenna  is  the  transducer  that  intercepts  propagating  electromagnetic  waves  and 
directs  them  into  guided  waves  in  a  transmission  line  In  the  case  of  a  transmitting 
antenna,  guided  electromagnetic  energy  in  a  transmission  line  is  converted  into  radiated 
electromagnetic  waves.  In  short,  the  antenna  couples  the  system  to  the  environment  [Ref 
1] 

Almost  all  radar,  communication,  and  electronic  warfare  (EW)  systems  require 
some  form  of  antenna  Even  for  highly  directive  antennas,  all  of  the  energy  is  not 
concentrated  into  the  main  beam;  there  is  always  radiation  in  directions  other  than  the 
intended  one  These  sidelobes,  as  they  are  called,  are  of  particular  importance  in 
electronic  warfare,  since  they  represent  a  major  limitation  to  the  system's  effectiveness 
Jamming  signals  or  false  targets  received  through  a  sidelobe  can  cause  confusion  and 
failure  of  a  mission  more  easily  than  false  signals  detected  by  the  main  beam.  The 
reduction  and  control  of  the  position  and  magnitude  of  sidelobes  is  extremely  important 
when  attempting  to  localize  target  azimuths  in  a  hostile  electronic  environment. 

Antenna  sidelobe  control  is  one  of  the  most  basic  methods  of  defeating  jammers 
and  avoiding  detection.  Low  sidelobes  also  reduce  multipath  signals,  which  can  cause 
signal  fading  in  communications  systems  and  tracking  errors  in  radars.  The  methods  of 


1 


sidelobe  suppression  for  passive  antennas  are  well  known,  and  the  technology  is  low  risk 
However,  low  sidelobes  are  accompanied  by  a  reduction  in  gain,  and  generally  require  a 
more  complicated  antenna  design  Consequently,  there  are  tradeoffs  in  the  antenna  design 
that  impact  the  overall  system  performance  and  cost 

B.  OBJECTIVE 

This  thesis  investigates  the  feasibility  of  reducing  antenna  sidelobes  for  the  Marine 
Radio  MRC-142.  This  system  provides  medium-range  ultra  high  frequency  (UHF) 
line-of-sight  multi-channel  radio  communications  over  a  distance  of  up  to  35  miles  The 
AN/MRC-142  radio  terminal  set  is  installed  in  the  M-998  variant  of  the  High  Mobility 
Multipurpose  Wheeled  Vehicle  (HMMWV)  and  when  deployed  appears  as  in  Figure  1 

This  unit  can  be  configured  as  a  repeater  to  extend  the  communications  range 
between  two  other  systems  or  to  overcome  terrain  conditions  that  prevent  a  direct 
line-of-sight  path  between  the  other  two  pieces  of  equipment  The  AN/MRC-142 
interfaces  with  a  distant  end  AN/MRC-142  over  a  1350  to  1850  MHz  UFIF  radio  link 
using  parabolic  antennas  that  are  mounted  on  masts  which  can  be  extended  up  to  50  feet 
The  current  design  uses  a  parabolic  reflector  with  a  focal-point  feed  and  has  a  first 
sidelobe  requirement  of  -10  dB  relative  to  the  mainbeam.  A  reduction  of  another  10  dB  is 
desired  for  the  new  design,  with  a  minimum  impact  on  the  system  upgrade  cost 
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50  Fool  Telescoping  Antenna 


M-998  H.MMWA 

Figure  1.  Deployed  AN/MRC-142  Radio  Terminal  Set  [Ref.  2]. 


C.  METHOD  OF  SOLUTION 

Possible  options  for  the  low  sidelobe  antenna  include 
I  an  array  antenna, 

2.  a  completely  redesigned  reflector  and  feed,  or, 

3  an  improved  feed  used  in  conjunction  with  the  current  reflector 
The  complete  antenna  specifications  are  given  in  the  Appendix  Options  I  and  2  could 
satisfy  the  specifications  in  their  entirety,  but  the  cost  of  design  and  development  would 
be  significant.  The  third  option  is  the  one  pursued  here  Even  though  all  of  the  antenna 
specifications  are  not  achieved  simultaneously  with  the  new  feed,  there  is  significant 
improvement  in  all  areas.  Most  importantly,  the  cost  of  the  new  feed  is  minimal  because 
of  its  simplicity. 
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The  proposed  feed  design  consists  of  a  fed  dipole  in  the  vicinity  of  a  parasitic  ring 
and  disk  as  shown  >n  Figure  2  [Ref  3]  Essentially  the  disk  acts  as  a  ground  plane  for  the 
dipole  and  *he  ring  narrows  the  beamwidth,  resulting  in  a  stronger  edge  taper  on  the 
reflector  The  dimensions  are  chosen  to  minimize  the  sidelobes  This  "tuning"  approach 
is  narrow  band  and,  as  discussed  in  Chapter  III,  does  not  allow  an  optimum  feed  design 
for  both  the  transmit  and  receive  frequencies 


Figure  2.  Dipole-Disk  Antenna  with  Ring  Used  as  a  Feed  in  Paraboidal  Reflector. 


D.  SCOPE  OF  THESIS 

In  Chapter  II  the  method  of  analysis  is  discussed  in  detail  All  of  the  gain  and 
pattern  data  presented  are  generated  using  the  method  of  moments  Chapter  III  discusses 
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the  tradeoffs  involved  in  the  antenna  design  As  mentioned  earlier,  the  greatest  challenge 
was  selecting  the  geometrical  parameters  to  provide  acceptable  performance  at  both 
frequencies  (1350  MHz  and  1850  MHz)  Finally,  Chapter  IV  presents  a  summary  of  the 
work  performed  along  with  some  conclusions 
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II.  THEORETICAL  BACKGROUND 


A.  E-FIELD  INTEGRAL  EQUATION  AND  THE  METHOD  OF  MOMENTS 

To  determine  the  pattern  and  gain  of  a  reflector  antenna,  it  is  first  necessary'  to  find 
the  current  on  the  reflector  surfaces  This  is  accomplished  by  solving  an  integral  or 
differential  equation  for  the  current  which  is  derived  from  Maxwell's  equations  and  the 
boundary  conditions  For  wire  and  surface  geometries  of  arbitral^'  complexity,  the 
integral  or  differential  equation  must  be  solved  numerically  on  a  high  speed  digital 
computer  The  method  of  moments  (MM)  is  such  a  technique  which  reduces  the  E-field 
integral  equation  (EFIE)  to  a  matrix  equation  which  is  then  solved  on  a  computer  [Ref  4] 
The  following  paragraphs  summarize  the  method  of  moments  technique  Phaser  notation 
is  used  throughout  and  e'"'  time  convention  is  assumed  and  suppressed 

If  the  current  on  a  body  is  known,  then  the  scattered  (or  radiated)  field  can  be 
determined  from  [Ref  5] 

E  =  -V(j)-j0)A  =  -j(0A-i;;^V(V«A),  (1) 

where  co  =  27rf  is  the  radian  frequency,  andp  and  e  are  the  permeability  and  peimittivity 
of  free  space  A  is  the  magnetic  vector  potential,  which  can  be  represented  as 

where  J ,  is  the  surface  current  and  S  refers  to  the  surface  of  the  conductor. 
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The  electric  field  integral  equation  is  based  on  the  boundary  condition  that  the  total 
tangential  electric  field  on  a  perfectly  electric  conducting  (PEC)  surface  of  an  antenna  or 
scatterer  is  zero.  This  can  be  expressed  as 

El(F=rs)  =  E{(r  =7s)  +  Ef(r  =rs)  =  0  onS  (3) 

where  S  is  the  conducting  surface  of  the  antenna  or  scatterer,  as  stated  before,  andr  =  r  s 
is  the  distance  from  the  origin  to  any  point  on  the  surface  of  the  antenna  or  scatterer  The 

subscript  "t"  indicates  tangential  components 

Referring  to  Figure  3,  Equations  (1)  and  (2)  can  be  combined  and  expressed  as 


Figure  3.  Coordinate  Systems  for  Computing  Radiation  Fields. 


(4) 


E^(?)  =  -jJlP'fUs(7')G(FJ')ds'  +  Vjf,V'.L(f')G(rJ')ds']- 

where 


and 


R=|?-r'|. 


(5) 


(6) 


Now,  applying  the  boundary  condition  of  Equation  (3)  gives 

P(?)  =  E;(r,)  =  f(P‘jjJ.(?')G(?,f')ds'  +  Vl|^V'.J,(f')G(r.r')ds']  (7) 

This  is  an  integral  equation  with  J  s  as  the  unknown  Equation  (7)  is  referred  to  as  the 
electric-field  integral  equation  In  Equation  (7).  V  and  V'  are  the  gradients  with  respect  to 

the  observation  (unprimed)  and  source  (primed)  coordinates,  respectively,  and 
G( r,r  ) 

is  the  Green's  function  for  a  three  dimensional  radiator  or  scatterer 

For  convenience,  an  operator  L  can  be  defined  which  will  operate  on  1. 

L(Js)=j(op|U,G,  ds'  +  '|v[v|JJ,G;  ds'  +||,J..VG;  ds'  (8) 

where  unknown  current  Js  depends  only  on  the  primed  quantities  G  is  a  scalar  quantity 
that  depends  on  both  the  primed  and  unprimed  quantities  It  is  desirable  to  have  all  the 

differential  operations  in  terms  of  the  primed  variables  of  integration  Equation  (8)  can 

therefore  be  written  as 


L(Js)=jcop||,J,G;  ds'  + 


V.J,ds'  +  jJ,J..V;  G;  ds'. 


(9) 
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where  the  vector  identity 


V(Js»  VG)  =  (Js*V)VG  +  (VG«  V)Js  +  Js  x(Vx  VG)  +  VGx(VxJ,)  (10) 

has  been  used 

The  gradient  of  G  in  terms  of  the  two  sets  of  coordinates  is  related  by 

VG  =  -0k  +  i)^r  =  -V'G  (11) 

Using  the  surface  divergence  theorem  developed  by  Mautz  and  Harrington  [Ref  5],  it  can 
be  shown  that 

11,  ^•(LG)ds'  =  0.  (12) 

Finally,  Equation  (7)  can  be  written  as 

M3.0  =  ll,ljcoTiIsG-l|V[(V'«Js)G]]ds'.  (13) 

Now  the  electric  field  integral  equation  can  be  represented  as 

E;(rs)  =  L(Js)  (14) 

Equation  (14)  can  be  solved  numerically  using  the  method  of  moments  The  current 
on  the  antennas  is  expanded  in  terms  of  a  finite  series  of  known  basis  functions  with 
unknown  amplitudes.  The  antenna  under  consideration  consists  of  surfaces  and  thin 
wires,  so  basis  functions  must  be  defined  on  each  On  a  surface  two  orthogonal 
components  are  required.  The  coordinate  system  for  the  reflector  is  shown  in  Figure  4 
and  the  following  basis  functions  defined  on  the  surface  [Ref  6]: 
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SEGMENTED  REFLECTOR  GENERATING  CURVE 


Figure  4.  A  Reflector  Geometry  for  Method  of  Moments. 

n  =  0,±l,...±x,  i=l,2 . Ns-2 

n  =  0,±I,...±x,  i=l,2 . Ns-1  '  ^ 

T,(t)  is  the  triangle  function  (which  extends  over  two  segments,  i  and  i  +  1)  and  P,(t) 
is  the  pulse  function  A  point  on  the  surface  of  the  antenna  is  specified  by  the  coordinates 
(t,  4)),  where  t  is  an  arclength  variable  along  the  reflector  generating  curve.  The  distance 
of  a  point  from  the  axis  of  symmetry  (z-axis)  is  given  by  p  and  N,  is  the  number  of 
surface  generating  points  A  discussion  leading  to  the  choice  of  these  basis  functions  is 
given  by  Mautz  and  Harrington  [Ref  5] 

For  thin  wires  the  radius  "a"  is  much  less  than  a  wavelength  and  the  current  will 

___  ^ 

only  have  an  axial  component.  Thus  the  operator  V  becomes  t(d/dt),  where  t  is  the  arc 
length  along  the  wire  Furthermore,  the  current  is  constant  around  the  wire  and  jj^ds' 


10 


reduces  to  2naiLdt',  as  shown  in  Figure  5  Thus  appropriate  basis  functions  on  the  wire 
are  [Ref  7]: 


=  i=l,2...,N..  (16) 

In  this  case  t  is  a  unit  vector  along  the  wire  axis  and  H  is  the  number  of  wire  generating 
points  The  total  current  is  the  sum  over  all  expansion  functions 


_  Nw-2  —  cc 

m=l  * 


N,-2  _ 

I  + 

p=,  ^=, 


I  lSnJ$n 


(17) 


The  coefficients  I  are  the  unknowns  which  need  to  be  determined 

From  the  linearity  of  L,  Equation  (14)  can  be  combined  with  Equation  (17)  as 
follows 


(18) 
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where  N^,  is  the  total  number  of  basis  functions  of  both  types  N^,  equations  need  to  be 

provided  to  enable  a  solution  for  the  unknown  constants  I,  Galerkin's  method  is  used  to 

reduce  Equation  (18)  to  a  matrix  equation,  and  weighting  functions  are  defined  such  that 
Wk  =  J  k  (Ref  8,  9]  Each  side  is  multiplied  by  each  of  the  basis  functions  and  then  each 
side  integrated  The  result  is  Nt,  equations  of  the  form 

E'.Wkds  =  [Jj,^,Wk*IiL(I,)  +  |l,^^Wk.l2L(J2)  +  ...+l^Wk*lNL(JN)]ds 
k=l,2,...N,  (19) 


which  can  be  written  as 


/f  E'«Wkds  =  /f  [Wk.Sl.L(J.)dsI 

■  k  •  ,_i 


(20) 


Interchanging  summation  and  integration. 


ff  E'.Wkds=  IlJlWk«L(J,)ds. 
■ «  1=1 


(21) 


Equation  (21 )  can  be  cast  in  the  form 


Vk  =  IZikl,. 


1=1 


(22) 


for  k  =  1,2,  where  the  impedance  elements  are  given  by 


Z,k  =  1 1,  ds  1 1  [jco pWk  •  J 1  -  i:(V '  •  J  .)(V  .  Wk)]Gds 


(23) 


The  excitation  elements  are  defined  by 


Vk  =  II,  Wk«E‘ds, 


(24) 
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In  matrix  notation  Equation  (22)  becomes 


[V1  =  [Z][I],  (25) 

and  the  unknown  coefficients  are  determined  by  solving  the  matrix  equation 

[I]  =  [Zr'[V]  (26) 

Submatrices  in  Equation  (25)  associated  with  the  t  and  $  surface  components  can  be 
identified  as  follows 


iv’n]  1 = r  [z;:][z‘„*]  ir  i‘  ■ 

[V*]  J  -  [  (zJ:*][zI:*]  JL  _  ■ 


(27) 


The  impedance  matrix  will  have  the  following  block  structure  when  all  of  the 
surface  and  wire  basis  functions  are  considered 


The  superscripts  s  and  w  refer  to  surface  and  wire,  respectively.  Although  n  in 
Equation  ( 1 7)  can  take  on  values  between  ±  ^c,  it  has  been  truncated  at  ±  M  in  the  above 
matrix.  In  practice,  when  the  feed  is  on  the  reflector  axis  of  symmetry,  a  converged 
solution  can  be  obtained  with  M  =  1.  If  the  pattern  is  calculated  for  the  antenna 


transmitting,  the  excitation  vector  will  have  only  one  non-zero  value  corresponding  to  the 
dipole  feed  point 

The  unknown  expansion  coefficients  are  determined  by  solving  the  matrix  in 
Equation  (27)  for  [I]  To  compute  the  electric  field,  the  coefficients  are  used  in  Equation 
( 1 7),  which  in  turn  is  used  in  Equations  ( 1 )  and  (2)  Define  a  measurement  vector, 

Rm  =jj,  Er«  Jsds,  (29) 

where  Er  denotes  a  unit  radiated  plane  wave  weighted  by  the  current  on  S  In  general  Is 
will  give  rise  to  both  0  and  4)  components  in  the  far  field  The  radiated  0  component  is  of 

the  form 

Er  =  =  0e-iMxu+yv+z«i  ^ 

where  u  =  sin0cos(t),  v  =  sin0sins(t),  and  w  =  cos0  By  substituting  Equation  (30)  into 
Equation  (29),  there  resjlts 

Rm=lj.se'J'‘'"‘'^''*^"'(0«Js)ds  .  (31) 

The  total  0  component  of  the  electric  field  is  the  superposition  of  all  surface  current 
contributions 

=  (32) 

The  analysis  holds  true  for  the  ({)  component.  The  measurement  elements  are 

Rm  =  jj,  js)ds  ,  (33) 
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and  the  electric  field  is  then  defined  by 

=  (34) 

^  4nr  m 

B.  GAIN  CALCULATION 

The  degree  to  which  any  particular  antenna  pattern  is  concentrated  into  a  beam  is 
known  as  directivity’  which  is  defined  by 

D  =  -^,  (35) 

where  is  the  maximum  power  flux  radiated,  and  P^vc  is  the  average  power  flux 
radiated  Since  the  average  power  flux  radiated  is  equal  to  the  total  power  radiated  (P^,,) 
divided  by  4rc  steradians, 

D  =  4rt^.  (36) 

Directivity,  as  defined  above,  is  a  function  of  the  antenna  pattern  only,  and  does  not  take 
into  account  the  antenna  efficiency  or  losses.  In  order  to  compare  the  relative 
performance  of  antennas  in  any  electronic-warfare  system,  it  is  necessary  to  consider 
antenna  efficiency.  The  pertinent  quantity  is  gain  G  which  is  defined  as 

G  =  KD,  (37) 


15 


where  K  is  the  efficiency  factor  (0^<  1)  Thus  it  is  seen  that  for  a  100%  efficient 
antenna  (K  =  1),  the  gain  would  be  equal  to  the  directivity  For  all  practical  antennas  K  < 
1  Gain  is  frequently  expressed  as  a  decibel  ratio 

G,db  =  101og,(jG.  (38) 

A  very  important  relationship  exists  between  gain,  physical  antenna  dimensions, 
and  electrical  wavelength,  expressed  by 

G  =  47t^,  (39) 

where  A  is  the  "effective"  area  of  the  antenna  and>.  is  the  wavelength  The  effective  area 
of  the  antenna  and  the  wavelength  must  be  expressed  in  the  same  units  (meters,  inches, 
feet,  etc  )  The  total  radiated  power  required  in  Equation  (36)  can  be  obtained  by 
integrating  the  electric  field  over  a  sphere 

p,„,  =  y|{|E„|^  +  |E*|^}r2sin0ded(t)  (40) 

0  0 

where  E,,  and  E,,,  are  known  from  Equations  (32)  and  (34) 
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III.  LOW  SIDELOBE  ANTENNA  DESIGN  STUDY 


A.  SUMMARY  OF  TRADEOFFS  IN  THE  ANTENNA  DESIGN 

The  purpose  of  this  thesis  is  to  achieve  a  reduction  in  the  antenna  sidelobes  for  the 
Marine  radio  AN/MRC-142  The  approach  will  be  to  upgrade  the  existing  reflector  in 
Figure  6  by  implementing  a  feed  design  similar  to  that  of  Kildal  and  Skyttemyr  [Ref  3] 
shown  in  Figure  2  of  Chapter  I 


Figure  6.  AN/MRC'142  Antenna  Assembly  and  Offset  Adapter  [Ref.  2.]. 


The  size  of  the  reflector  aperture  is  36  x  48  inches  as  shown  in  Figure  7  The 
projected  aperture  is  not  circular,  but  has  two  circular  segments  removed.  The  method  of 
moments  code  that  was  used  to  compute  the  gain  and  radiation  pattern  can  only  handle 
rotationally  symmetric  reflectors  In  the  design  of  the  antenna,  circular  parabolic 
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reflectors  of  36  and  48  inches  were  analyzed  Additionally,  the  structure  supports  only 
vertical  currents,  whereas  the  method  of  moments  code  assumes  radial  and  azimuthal 
currents  (solid  surface)  The  actual  AN/MRC-142  performance  will  lie  somewhere 
between  the  two  The  pattern  will  resemble  the  48  inch  reflector  in  the  plane  of  the  wide 
dimension,  and  the  36  inch  reflector  in  the  plane  of  the  narrow  dimension  The  effect  on 
gain  can  be  accurately  estimated  based  on  the  area  reduction  from  the  48  inch  diameter 
circular  antenna  due  to  removal  of  the  circular  segments 


Figure  7.  Antenna  Radius  Determination. 


The  reflector  surface  is  a  section  of  a  parabaloid  and  is  fed  from  the  front  Four 

factors  affect  the  efficiency  of  this  type  of  antenna: 

1  Spillover  efficiency, 

2  Illumination  efficiency, 

3  Blockage  by  the  feed; 
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4  Cable  and  line  losses  from  the  antenna  input  to  the  feed 


Taking  into  consideration  all  these  factors,  most  of  the  burden  in  maximizing  the 
efficiency  is  placed  upon  the  antenna  feed  The  dipole-disk  antenna,  which  essentially 
consists  of  a  halfwave  dipole  over  a  circular  plane  reflector,  is  w  idely  used  as  feed  for 
paraboloidal  reflector  antennas  at  MiF  and  UHF  A  method  of  improving  this  feed 
consists  in  incorporating  a  circular  conducting  beam-forming  ring  (BFR)  over  the  dipole 
in  a  plane  parallel  to  the  disk  The  beam-forming  ring  compresses  the  H-plane  pattern  of 
the  dipole-disk  feed,  whereas  it  has  no  significant  effect  on  the  E-plane  pattern  U  hen 
this  antenna  is  used  as  feed  for  paraboloidal  reflectors,  it  exhibits  near-identical  principal 
plane  aperture  illumination  and  reduces  H-plane  spillover  lobes 

The  diameter  of  the  disk  should  be  as  small  as  possible  in  order  to  reduce  blockage 
loss  However  a  small  disk  and  ring  result  is  a  broad  feed  radiation  pattern  which  does 
not  provide  sufficient  edge  taper  for  a  low  sidelobe  pattern  To  obtain  the  best 
performance  a  parametric  study  was  performed  for  the  geometrical  quantities  shown  m 
Figure  8  and  listed  in  Table  1  The  parameters  and  their  effect  in  antenna  radiation 
pattern  and  gain  discussed  in  the  following  sections  All  dimensions  used  by  the 
computer  code  are  in  wavelengths 
I.  Reflector  Size 

The  size  of  the  reflector  is  36  x  48  inches  as  shown  in  Figure  7  Two  different 
designs  were  obtained  for  each  frequency:  one  for  36  inches  and  the  other  one  for  48 
inches  At  1350  MHz,  the  wavelength  is  8  75  inches  and  at  1850  MFIz  the  wavelegth  is 
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6  38  inches  Thus  36  inches  corresponds  to  4  IX  and  5  6X,  while  48  inches  corresponds 
to  5  5X  and  7.5X  at  these  two  frequencies 

The  focal-length-to-diameter  ratio  (f/D)  of  the  existing  reflector  antenna  could 
not  be  determined  from  AN/MRC-142  technical  manuals  Typical  values  are  in  the  range 
of  0.3  to  0.55,  Smaller  values  are  preferable  to  reduce  spillover  loss  and  increase  edge 
taper  A  value  of  0  4  was  chosen  for  f/D  This  parameter  could  be  varied  slightly  without 
significant  modification  of  the  antenna's  gain  and  pattern 


Figure  8.  Program  Parameters 
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DM  = 
DS  = 
ZC  = 
FOD  = 
ZR  = 
RC  = 
SW  = 
DLN  = 
PLN  = 
DDIP  = 
AW  = 
HLGTH  = 
DHP  = 


TABLE  1.  DESIGN  PARAMETERS _ 

diameter  of  paraboloid 
diameter  of  disk 

distance  of  dipole  to  ground  plane 
focal  length  to  diameter  ratio  of  the  parabolid 
distance  of  ring  from  fed  dipole  along  the  z  axis 
radius  of  ring 

ring  width  (modeled  as  a  thin  strip) 

fed  dipole  length 

parasitic  director  dipole  length 

distance  between  DLN  and  PLN  along  the  z  axis 

radius  of  dipoles 

parasitic  (H  plane)  dipole  length  (not  used) 

H  plane  dipole  distance  from  feed  along  the  y  axis  (not  used) 


2.  Gain  Calculation 

Gain  estimates  for  the  reflector  and  feed  were  obtained  by  first  integrating  the 
total  electric  field  to  find  the  total  radiated  power  The  result  is  then  used  in  Equation 

[36]  to  determine  the  directivities  shown  in  Table  2.  These  numbers  include  spillover  and 
aperture  efficiency  because  the  total  field  is  integrated  (that  is,  feed  plus  scattered  fields). 

Losses  that  are  not  accounted  for  and  thus  determine  the  value  of  K  in  Equation 

[37]  include: 

1  Transmission  line  losses  for  the  cable  that  feeds  the  dipole.  Assuming  0.5  dB/ft 
and  3  ft  of  cable  yields  1 .5  dB  of  loss. 
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2  Reflection  loss  from  the  reflector  surface  This  is  due  to  the  fact  that  the  surface  is 
a  grid  rather  than  a  solid  A  conservative  estimate  is  0  5  dB 


3  Any  impedence  mismatch  loss  A  VSWR  of  2:1  results  is  approximately  0  5  dB 
of  loss 


Based  on  the  above  numbers  the  total  loss  is  about  2  5  dB  The  corresponding 
gains  are  shown  in  Table  3 


TABLE  2 

.  CALCULATED  DIRECTIVITIES. 

48  inches 

36  inches 

1350  MHz 

24.75  dB 

22.19  dB 

1850  MHz 

27.44  dB 

24.90  dB 

TABLE  3. 

CALCULATED  GAINS  WITH  LOSS. 

48  inches 

36  inches 

1350  MHz 

22.25  dB 

19.19  dB 

1850  MHz 

23.94  dB 

21.90  dB 

B.  MODIFICATIONS  TO  ANTENNA  DESIGN 

Each  of  the  parameters  of  the  feed  and  reflector  were  varied  in  small  increments  and 
their  effect  on  the  antenna  pattern  and  gain  was  evaluated.  The  results  are  summarized 


below. 


1.  Baseline  Feed  Design 

The  dimensions  of  the  baseline  feed  design  are  listed  in  Table  4  The  disk  serv  es 
as  a  ground  plane  for  the  dipole.  The  phase  center  of  the  feed  will  vary  with  frequency 
and  angle  If  the  disk  is  an  effective  ground  plane  the  phase  center  of  the  dipole  and  its 
image  will  be  approximately  at  the  midpoint  of  the  two.  However  in  the  following 
discussion  the  dipole  is  considered  to  be  the  center  of  the  feed 

Low  sidelobes  require  a  strong  amplitude  taper  across  the  reflector,  which  in  turn 
requires  a  narrow  feed  beamwidth  Generally  a  narrow  feed  beamwidth  is  associated  with 
a  feed  aperture  which  results  in  increased  blockage  Thus  there  is  a  tradeoff  between 
sidelobe  taper  and  blockage  loss 

The  geometrical  parameters  in  Table  4  were  the  starting  point  for  the  parametric 
study,  which  is  summarized  in  the  following  sections.  The  baseline  antenna  performance 
is  shown  in  Figure  10.  The  solid  line  is  the  E-plane  (4)=0°)  pattern  and  the  dashed  line  the 
H-plane  ((})=90°)  pattern. 
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Relative  Power,  dB 


Figure  10.  Baseline  Antenna  Pattern. 


TABLE  4.  BASELINE  FEED  GEOMETRY 

Z  distance  or 

Variable  Value  coordinate’ 

disk  diameter 

4.1 

Zc  =  .25 

ring  diameter 

0.7 

Zr  =  0.0 

fed  dipole  length 

0.5 

Zd  =  -0.25 

parasitic  dipole  length 

not  part  of  baseline  design 

fed  dipcie  height  above  disk 

0.5 

Zc-Zd 

ring  heigth  above  disk 

0.25 

Zc-Zr 

wire  radius 

.OOU 

- 

'  The  fed  dipole  is  located  at  Z  =  0;  positive  Z  is  away  from  the  reflector 
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2.  Disk  Size  and  Location 


The  disk  must  be  large  enough  to  serve  as  an  effective  ground  plane  for  the 
dipole.  Generally  the  dipole  height  must  be  significantly  less  than  the  diameter  of  the 
disk  For  the  dipole  and  its  image  to  add  constructively  in  the  direction  of  the  reflector, 
the  spacing  should  be  approximately  25X  The  patterns  for  disk  locations  of  Z  =  +  IX, 
and  -  IX  are  shown  in  Figures  11  and  12  All  other  feed  parameters  are  held  constant 
(values  given  in  Table  4)  The  reflector  diameter  is  4  IX  and  f/D  =  05  in  all  cases.  The 
high  levels  in  the  H  plane  patterns  for  80°  <  0  <  120°  are  due  to  feed  spillover  A 
comparison  of  the  figures  shows  that  there  is  a  tradeoff  between  spillover  levels  in  the 
rear  hemisphere  (90  <  0  <  180)  and  sidelobe  levels  in  the  forward  hemisphere  {0<  90°) 
The  impact  of  changing  the  disk  radius  by  .  IX  is  illustrated  in  Figures  13  and  14 
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Relative  Power,  dB 


Figure  11.  Radiation  Pattern  with  Disk  Shift  of  +0.1  X. 
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Relative  Power,  dB 


Figure  12.  Radiation  Pattern  with  Disk  Shift  of  -0.1  X. 
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Relative  Power,  dB 


Figure  13.  Radiation  Pattern  with  Disk  Size  .iX. 
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Relative  Power,  dB 


Figure  14.  Radiation  Pattern  with  Disk  Reduced  -.1X. 


29 


3.  Ring  Size  and  Position 

The  beam-forming  ring  (BFR)  compresses  the  H-plane  beamwidth  but  leaves  the 
E-plane  beamwidth  essentially  unaffected  Thus  the  ring  size  and  location  can  be  used  to 
equalize  the  feed  beamwidth  Patterns  for  two  ring  radii  are  shown  in  Figures  15  and  16 
The  ring  acts  somewhat  like  a  director  in  a  Yagi  array  when  it  is  located  in  front 
of  the  dipole  (towards  the  reflector)  [Ref  6]  When  the  ring  is  located  behind  the  dipole  it 
has  little  effect  on  the  reflector  scattered  field  As  illustrated  in  Figure  17  and  18,  a  shift 
of  R  results  in  about  a  1  dB  difference  in  the  first  sidelobe  level 
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Relative 


Figure  15.  Radiation  Pattern  with  Ring  Size  of  0.6X.. 
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Figure  16.  Radiation  Pattern  with  Ring  Size  of  O.SX. 
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Figure  17.  Radiation  Pattern  with  Ring  at  Z  =  -0.2X. 
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0 


Figure  18.  Radiation  Pattern  with  Ring  at  Z  =  -0.1  X. 

4.  Adding  a  Parasitic  Dipole 

A  parasitic  dipole  was  added  as  a  director  in  order  to  improve  the  directivity  of 
the  feed  and  therefore  increase  the  reflector  edge  taper  This  method  is  effective  in 
narrowing  the  beamwidth  of  cavity-backed  dipole  antennas  In  this  case  it  was  found  that 
changing  the  size  and  position  of  the  parasitic  element  has  no  significant  effect  on  the 
performance  of  the  antenna  and  does  not  offer  any  advantage.  This  is  probably  because 
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the  parasitic  dipole  is  not  located  within  a  cavity  Adding  a  second  ring  around  the 
parasitic  dipole  could  improve  the  performance  In  an  effort  to  keep  the  feed  as  simple  as 
possible,  the  parasitic  element  was  not  considered  further 

5.  Defocusing  the  Feed 

The  phase  center  of  the  feed  varies  with  frequency  and  angle,  but  is  generally 
located  between  the  fed  dipole  and  disk  The  antenna  bandwidth  can  be  increased  by 
slightly  displacing  the  feed  from  the  reflector  focus  Defocusing  can  cause  a  degradation 
at  one  frequency  but  an  improvement  at  the  other  Relative  to  an  optimum  pattern, 
defocusing  will  result  in  a  spoiled  beam;  that  is,  null  filling  and  change  in  beam  slope  An 
example  is  shown  in  Figure  19  Defocusing  also  affects  pattern  levels  in  the  rear 
hemisphere. 

C  OPTIMUM  DESIGNS 

I.  Single  Frequency  (1850  MHz) 

After  making  all  the  comparisons  mentioned  earlier,  an  optimum  configuration 
was  identified  for  the  36  inch  parabolic  reflector  operating  at  1850  MHz  The  patterns  are 
shown  in  Figure  20,  and  the  configuration  parameters  shown  in  Table  5 
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Relative  Power,  dB 


Relative  Power,  dB 


Figure  20.  Radiation  Pattern  of  Optimum  Single  Frequency  Design  for  36  inch  Reflector 
(The  Dimensions  are  given  in  Table  6).  _ 


This  radiation  pattern  meets  almost  all  the  requirements  listed  in  the  Appendix 
A  The  calculated  total  gain  is  G  =  20  7775,  which  is  sufficient,  but  the  far  out  sidelobes 
are  about  2  dB  too  high  Unfortunately  this  is  a  narrow  band  design  and  has  poor 
performance  at  the  lower  frequency  (1350  MHz)  This  design  illustrates  that  achieving 
the  specifications  at  a  single  frequency  is  relatively  easy.  However,  broadbanding  the 
antenna  proves  to  be  much  more  difficult 

2.  Dual  Frequency  (1350  and  1850  MHz) 

Starting  from  the  configuration  described  in  Table  4,  attempts  have  been  made  to 
obtain  a  broadband  design  for  the  frequency  range  from  1350  to  1850  MHz  Figures  21 
and  22  present  the  radiation  patterns  obtained  for  the  36  inch  antenna  which  achieves 
most  of  the  specifications  given  in  the  Appendix  A.  The  gains  are  marginal  at  both 
frequencies  and  the  beamwidth  at  the  low  frequency  in  too  wide.  Table  6  describes  the 
configuration  of  this  reflector. 
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Figure  21.  Radiation  Pattern  of  36  inch  Reflector  for  Lower  Frequency  (1350  MHz). 
Dimensions  are  given  in  Table  6. _ 
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Relative  Power,  dB 


Figure  22.  Radiation  Pattern  of  36  inch  Reflector  for  Higher  Frequency  (1850  MHz). 
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TABLE  6.  OPTIMUM  GEOMETRY  FOR  36  INCH 
REFLECTOR  FOR  THE  LOW  FREQUENCY 

DM 

a 

4.1 

DNL 

a 

0.5 

DS 

a 

0.6 

PLN 

a 

0.4 

ZC 

a 

0.3 

DDIP 

a 

-0.25 

FOD 

a 

0.4 

ROT 

a 

0 

RC 

a 

0.6 

AW 

a 

0 

ZR 

a 

0.1 

ZSHIFT 

a 

0.25 

SW 

a 

0.05 

For  the  48  inch  antenna.  Figures  23  and  24  show  the  radiation  patterns  for  the 
optimum  configuration  whose  dimensions  are  listed  in  Table  7 
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Figure  23.  Radiation  Pattern  of  48  inch  Reflector  for  Lower  Frequency  (1350  MHz). 
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Figure  24.  Radiation  Pattern  of  48  inch  Reflector  for  Higher  Frequency  (1850  MHz). 
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TABLE  7.  OPTIMUM  GEOMETRY  FOR  48  INCH 

REFLECTOR  FOR  HIGH  FREQUENCY 

DM 

an 

DNL 

0.48 

DS 

a 

0.8 

PLN 

0.4 

ZC 

a 

0.3 

DDIP 

-0.25 

FOD 

a 

0.4 

ROT 

0 

RC 

a 

0.5 

AW 

0.001 

ZR 

a 

0 

ZSHIFT 

0.2 

SW 

a 

0.05 

It  is  worth  noting  that  there  is  no  feed  configuration  that  can  be  used  for  both  frequencies 
and  simultaneously  meet  the  specifications 
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rv.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  objective  of  this  thesis  was  to  find  the  optimum  combination  of  the  feed  and 
reflector  geometry  over  the  frequency  band  1350  to  1850  MHz  The  current  antenna  size 
is  36  X  48  inches  Since  the  computer  code  only  accepts  circular  reflector  apertures,  two 
separate  circular  antennas  were  analyzed  One  was  36  inches  and  the  second  was  48 
inches.  The  proposed  feed  design,  which  consists  of  a  fed  dipole  in  the  vicinity  of  a 
parasitic  ring  and  disk,  was  a  good  approach  to  achieve  a  significant  improvement 
compared  with  the  original  geometry  of  the  Marine  radio  AN/MRC-142.  Because  of  its 
simplicity,  the  cost  of  the  new  feed  represents  a  minimum  impact  on  the  system  upgrade 

Calculations  for  both  sizes  were  made  separately  and,  as  it  is  shown  in  Chapter  III, 
the  36  inch  reflector  could  not  satisfy  the  requirements  at  both  frequencies 
simultaneously,  but  it  gave  the  best  result  for  the  high  frequency  (1850  MHz)  The  results 
obtained  with  the  48  inch  reflector  come  very  close  to  satisfying  all  the  requirements, 
although  the  sidelobe  level  (SLL)  is  off  by  approximately  2  dB  at  wide  angles  in  the 
H-plane 

An  assumption  with  regard  to  the  transmission  line  losses  was  made,  which  affect 
the  total  gain  of  the  antenna.  In  the  field  calculations,  solid  surface  reflection 
characteristics  were  assumed.  In  fact,  the  surface  is  a  grid,  and  consequently  the 
reflection  coefficients  for  the  0  and  (j)  components  of  the  electric  field  will  vary  with  the 
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angle  of  incidence  This  effect  will  be  small  if  the  wire  spacing  is  small  compared  to  X. 
Achieving  the  required  gain  was  not  a  problem  for  the  48  inch  design,  and  small 
adjustments  in  the  antenna  feed  geometry  did  not  have  a  significant  negative  impact  on 
the  antenna  performance.  It  was  also  determined  that  if  there  was  an  improvement  due  to 
an  additional  parasitic  dipole,  the  increase  in  complexity  did  not  merit  its  use 

To  summarize,  a  low  sidelobe  reflector  antenna  has  been  designed  which  fulfilled 
most  of  the  requirements  for  the  Marine  radio  The  dipole  disk  feed  with  beam-forming 
ring  is  electrically  and  mechanically  simple  It  can  be  supported  by  means  of  a  stiff 
coaxial  cable  aligned  along  the  axis  of  the  paraboloid,  and  held  by  a  center  support  which 
offers  minimum  blockage  when  compared  with  the  usual  strut  supports,  as  was  stated  by 
Kildal  and  Skyttemyer  [Ref  3] 
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APPENDIX  A.  ANTENNA  SPECIFICATIONS 


Over  the  frequency  range  of  1350-1850  MHz,  the  antenna  must  have  gains  as 
described  in  Figure  A. 


Azimuth  Relative  to  Boresight  (degrees) 

Gain  (dBi) 

AtO 

A  -  see  note  A 

At±  10 

At  least  5  dB  lower  than  A 

At  ±20 

At  least  19.5  dB  lower  than  A 

At  ±30  to  180 

At  least  27  dB  lower  than  A 

Shall  satisfy  the  relationship  17.5  <  A  <  21  dB 


Figure  A.  Gain  Requirements. 
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Frequency  Range 


Gain 


TABLE  A.  ANTENNA/MAST  REFERENCE  DATA 


{Lightweight,  open  grid,  reflector 


1350  to  1850  MHz 


1350  MHz 


1850  MHz 


Breamwidth 


Horizontal 


IV2  ±  1  °  from  boresight 


1 5  ±  1  °  from  boresight 


Sidelobes 


Antennna 


Vertical 


7’/?  a:  1°,  -2°  from  boresight 


15  ±  1°,  -2°  from  boresight 


Sidelobes 


Polarization 


Weigth 


Maximum  Input  Power 


Maximum  Output  Power 


Reflector  Dimensions 


Mast 


Height 


Length  (for  storage  and  transit) 


Maximum  Height  (deployed) 


Headload 


Maximum  Wind  (50  ft.  erection) 


Weight 


2-Man  Deployment  Time 


Mast  Accessories  Time 


18  dBi  minimum 


20  dBi  minimum 


-h7.5° 


3  dB  half  power 


15  dB  down 


10  dB  min 


3  dB  half  power 


15  dB  down 


10  dB  min 


Horizontal  or  vertical  polarization 


17.5  lbs 


15  W 


2320  W  Effective  Radiated  Power  (ERP) 


45  in  X  35  in.  Depth;  16  in 


Telescoping  mast,  consists  of  eleven 
carbon-fiber  tubular  sections 


Continously  adjustable  from  5  ft  8  in  to  50  ft 


5  ft  8  in 


50  ft 


37  lbs  maximum 


60  mph  (without  ice) 

40  mph  (with  Vz  in  ice  loading) 


35  lbs 


15  minutes 


Stowage  of  mast  accessories 


Total  accessory  bag  weight:  44  lbs 
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APPENDIX  B.  PROGRAM  CODE* 


QamZf 

C  PROGRAM  mrcint.f 
C 

C  FAR  FIELD  PATTERN  INTEGRATION  FOR  GAIN 
CALCULATIONS. 

C  ganaral  wira/bor  gaomatrias  can  ba  handlad  *'* 

C 

C  »»  INPUT  FILE  IS  pcurranl  -  SPECIALIZED  FOR 
PROGRAM  mrc.f 

C  »»  HAS  A  GAIN  NORMALIZATION  LOOP 
INCORPORATED 

C  »»  WIRE  DOES  NOT  HAVE  TO  LIE  IN  THE  Z.0 

PLANE 

C 

DIMENSION 

RH(1000)7H(1000),AG(100),XG(1CIO),XT(10),AT(10) 

DIMENSION 

Q(30).S<30).nw1(2),nw2(2),pha(1000).xh(1000).yh(1000) 

COMPLEX 

RS(1800),ETW,EPW,U,CNT(1500).ET1  ,EP1  .ET2,EP2 
COMPLEX 

CONJG,CMPLX.EXP1,EXP2.RW(2000),ET.EP,UC 
DATA  pi>3.141S9/ 
data  nl/2/.xt(1),af(1)/.S773S03,1./ 
c  data 

nt'4/,i(t(  1 )  .xt(  1 0)  .at(  1 )  ,at(  1 0)/.33998 1 04 .  .86 1 1 363 1 1 5 , 
c  '  .6S21451S48..34785484S1/ 

RAD=PI/180. 

BK=2.'PI 

U.(0..1.) 

UC«-U/4./PI 

OPEN(1.FILE«outgaus') 

READ(1.*)  NPHI 
NT2=NT/2 
D01  K»1.NT2 
K1=NT-K+1 
AT(K1)=AT(K) 

XT(K1)=XT(K) 

1  XT(K)»-XT(K) 

D0  3K=1,NPHI 
3READ(1.*)  XG(K).AG(K) 

OPEN(3.fila«'pcurrent', STATUS*  old') 

READ(3.-) 

DM.AWZC.RC.DS.NP1,NP2,NBLOCK.NWIRES.FFAC 

READ(3.*) 

NP.MT1  .MP1  .N,MT2.NSURF,NROW.Mrll, MODES 

writa(6.*)  ’NP1  ,NP2,N,MODES.NROW,NWIRES,NP=' 

writa(6,")  NP1  ,NP2,N, MODES, NROW.NWIRES,NP 

READ(3.*)  (xh(i),i-1.rp) 

READ(3,*)  (yh(i),i-1,np) 

READ(3,‘)  (zh(i),i=1,np) 

READ(3,')  (rti(i),i«1,np) 

READ(3,‘)  (pha(i),i=1,np) 

READ(3,*)  (nw1(i).i»1,nwira8) 

READ(3,*)  (nw2(i),t»1,r!wiras) 


READ(3.*)  (cnl(i).Kl.rwow) 

CLOSE(3) 

NS*NBLOCK‘N 

OPEN(8.FILE>'oiitint') 

WRrrE(8.8000)  NP1.NP2.NBLOCK 
8000  FORMAT!//, 10X.’**  WIRE/BOR  ANTENNA  GAIN 
CALCULATION 

•//,10X,'NUMBER  OF  BOR  POINTS  (NP1).'.i4,/,10x. 
‘■NUMBER  OF  DIPOLE  POINTS  (NP2)«',I4,/.10X. 
‘■NUMBER  OF  MODE  BLOCKS  (2’M+1)>'.I2) 

NDIVT-2 

NDIVP-2 

WRrTE(6.‘)  ENTER  NDIVT  AND  NDIVP^ 

READ(5.‘)  NDIVT,NDIVP 
S1>0. 

S2>180 

01=0. 

02*180. 

vmrrE(6,‘)  theta  start  and  stop  angles 

(DEGREES)' 

READ(5.‘)  S1.S2 

WRITE(6.‘)  THETA  START  AND  STOP  ANGLES 
(DEGREES)' 

READ(5.‘)  01.02 
C 

C  DEFINE  STANDARD  BOR  INDICES 
C 

MT1=NP1-2 

MP1.NP1-1 

MT2.NP2-2 

C  NWIRES=1  SINCE  DUMMY  SEGS  WILL  HAVE  ZERO 
CURRENT 
NWIRESsI 
NW1(1)«NP1+1 
NW2(1)=NP 
C 

C  GAIN  NORMALIZATION  LOOP  (USE  THETA«PHI=0) 

C 

ET1=(0..0.) 

ET2»(0..0.) 

ETW=(0..0.) 

CALL 

PLANEW(NWIRES.NW1  ,NW2,NP1  ,NP2,xh,yh,zh.0..0.,RW) 
DO610L-1.MT2 

610  ETW=ETW+RW(L)‘CNT(L+NS) 

EXP1=(1.,0.) 

EXP2»CONJG(EXP1) 

DO  611  M>1.MHI 
NM=M-1 

CALL  PLANES(NM.NP1 .0,NT,RH7H,XT,AT,0.,RS) 
NTOP1-MODES-NM 
NTOP2»NBLOCK-(NTOP1 +1 ) 

NS2»NTOP1‘N 
NS1*NTOP2‘N 
D0  612L»1,MT1 

ET1»ETURS(L)‘CNT(L*NS1  )‘EXP1 
IF(NM.EO.O)  GO  TO  612 
ET1.ETURS(L)‘CNT(L*NS2)‘EXP2 


*  Although  this  is  FORTRAN  code,  the  column  conventions  are  not  observed  for 
this  listing  so  that  space  can  be  conserved. 
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612  CONTINUE 
00  613  L>1.MP1 

ET2«ET2+RS(UMT1)‘CNT(L+NS1+MT1)*EXP1 
IF(NM.EQ.O)  GO  TO  613 
ET2»ET2-RS(L+MT1  )‘CNT(L*NS2+Kin'1  )-EXP2 

613  CONTINUE 
611  CONTINUE 

ET.UC*(ET1  ♦ET2+ETW) 

EMAX>CABS(ET) 

VVRITE(6,')  'EMAX.'.EMAX 
WRITE(8,8008)  NT.NPHI.NDIVP.NDIVT.SI  .32.01 .02 
8008FORMAT(//.10X,'NT«M3./.10x.'NPHI=M4./.10X. 

•  NUMBER  OF  INTEGRATION  INTERVALS  IN 
PHI='.i3,/.10x. 

•  NUMBER  OF  INTEGRATION  INTERVALS  IN 
THETA=’,i3,/.10x. 

•  THETA  LIMITS.  START,STOP='.F7.2,2X.F7.2./.10X. 

•  'PHI  LIMITS:  START.STOP.',F7.2.2X,F7.2) 
c 

c  thata  intagratton  intarvate 
c 

S1-SrRAD 

S2=S2‘RAD 

Q1=Q1*RAD 

Q2=Q2'RAD 

DS=(S2-S1  )/FLOAT(NDIVT) 

DO  502  l=1.NDIVT+1 
S(I)»=(I-1)-DS 
WRrrE(6.-)  •|.S(I)=M,S(I) 

502  CONTINUE 
c 

c  phi  intagration  inlarvals 
c 

DQ.(Q2-Q1)/FLOAT(NDIVP) 

00  503  U1.NDIVP+1 
Q(l)«(l-ir0Q 
WRITE(6.*)  •|.Q(I).M.Q(I) 

503  CONTINUE 
SUM=0. 

do2000JJ=1.NDIVP 

P1=DQ.'2. 

P2=(Q(JJ+1)+Q(JJ))/2. 

00  2000  J=1.NPHI 
PHR=prxg(j)+p2 
DO  1000  IU1.NDIVT 
T1=:DS/2. 

T2=(S(ll+1)+S(ll))/2. 

DO  1000  Isl.NPHI 
THR=T1*XG(I)+T2 
IFfTHR.LT.PI)  GO  TO  510 
THR=(BK-THR) 

PHR=PHR+PI 
510  CONTINUE 
STRA=SIN(THR)*AG(l)’ag(j) 

ET1=(0.,0.) 

EP1=(0..0.) 

ET2=(0..0.) 

EP2»(0..0.) 

ETW=(0,.0.) 

EPW=(0.,0.) 

C 

C  WIRE  CONTRIBUTION  TO  THE  RADIATED  FIELD 


CALL 

PLANEW!NWIRES.NW1  .NW2.NP1  .NP2.XH.VH2H.THR,P 
HR.RW) 

DO210L«1.MT2 

ETW.ETW*RW(L)-CNr(L+NS) 

210  EPW»EPW+RW(L+MT2)‘CNT(L+NS) 
c 

c  BOR  SURFACE  CONTRIBUTION  TO  THE  RADIATED 

FIELD 

c 

D0  300M.1.MHI 
NM=M-1 

EXP1»CEXP(CMPU(O..FLOAT(NM)'PHR)) 

EXP2«CONJG(EXP1) 

CALL  PLANES(NM.NP1 .0.NT.RH2H.XT.AT.THR.RS) 
NTOPUMODES-NM 
NTOP2«NBLOCK-(NTOP1+1 ) 

NS2«NTOPrN 

NS1=NTOP2'N 

DO250L-1.MT1 

ET1«ETURS(L)'CNT(L+NS1  )*EXP1 
EP1-EP1+RS(L+N)’CNT(L+NS1)*EXP1 
IF(NM.EQ.O)  GO  TO  250 
ET1.ET1+RS(L)-CNT(L+NS2)‘EXP2 
EP1=EP1-RS(L+N)*CNT(L+NS2)‘EXP2 
250  CONTINUE 
DO260L=1.MP1 

ET2=ET2+RS(L+MT1  )’CNT(L+NS1  +MT1  j’EXPI 
EP2=EP2+RS(L+MT1+N)‘CNT(L+NS1+MT1)‘EXP1 
IF(NM.EQ.O)  GO  TO  260 
ET2=ET2-RS(L+MT1)*CNT(L+NS2+MT1)*EXP2 
EP2=EP2+RS(L+MT1  ♦N)‘CNT(L+NS2+MT1  )-EXP2 
260  CONTINUE 
300  CONTINUE 
ET»UC*(ET1+ET2+ETW) 

EPaUC'(EPUEP2+EPW) 
SUM.SUM4.(CABS(ET)**2+CABS(EP)”2)*STRA 
1000  CONTINUE 
2000  CONTINUE 
PRAD=TrPrSUM 
GAIN=4.‘PrEMAX"2T>RAD 
CDB=10.'alog10(GAIN) 

VI/RnE(8.8010)  PRAD.EMAX.GAIN.GDB 
8010  FORMAT(//.10X.TOTAL  RADIATED 
POWER»'.E15.8./.10X. 

•  -MAXIMUM  FIELD  VALUE  (V/M)=’.E15.8./.10X. 
••GAIN='.E15.8./.10x.'(dB)='.l8.2) 

WRnE(6.’)  'PRAD.GAIN.dBs’.PRAD.GAIN.GDB 
9999  STOP 
END 

SUBROUTINE 

PLANES(M1  .NP1  ,NP2.NT.RH2H.XT.AT.THR.R) 

C 

C  PLANE  WAVE  EXCITATION  VECTOR  FOR  BOR 
ELEMENTS. 

C  NO  CHANGE  FROM  HARRINGTON'S  ORIGINAL 
PROGRAM  OTHER  THAN 
C  ••  ONLY  ONE  MODE  PER  CALL  " 

C 

COMPLEX 

R(1800).U.U1.UA.UB.FA(50).FB(50).F2A.F2B.F1A.F1B 
COMPLEX  U2.U3.U4.U5,CMPLX 
DIMENSION  RH(  1 000)2H(  1 000), XT(  1 0).AT(  10) 
DIMENSION  R2(10)Z2(10).BJ(5000) 

U.(0..1.) 
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U1«3.141593*U”M1 

NF-1 

M2-M1 

NP-NP1+NP2 

MP-NP-1 

MT.MP-1 

N>MT+MP 

N2.2-N 

CC»COS(THR) 

SS«SIN(THR) 

M3>M1+1 

M4>U2+3 

IF(MI.EQ.O)  M3«2 

MS«MU2 

M6>M2'f2 

DO  12  IP-1, MP 

K2-IP 

I.IP-»1 

DR-.5’(RH(I)-RH(IP)) 

DZ-.S‘(ZH(I)-ZH(IP)) 

01=SQRT(DR-DR+DZ'DZ) 

R1«.25-(RH(I)+RH(IP)) 

IF(ABS(R1).LT.1.E-5)R1=1. 

Z1=.S‘(ZH{I)+ZH(IP)) 

DR=.5"DR 
D2-DR/R1 
DO  13L-1.NT 
R2(L)xR1+DR"XT(L) 
Z2(L)=ZUDZ'XT(L) 

13  CONTINUE 
D3=DR'CC 
D4=-DZ‘SS 
D5=DrCC 
DO  23  MxM3.M4 
FA(M)-0. 

FB(M).0. 

23  CONTINUE 
DO  1SL-1,NT 
X-SS'R2(L) 

IF(X.GT..5E-7)  GOTO  19 
DO  20  MxM3.M4 
BJ(M)=0. 

20  CONTINUE 
BJ(2)=1. 

S=1. 

GOTO  18 
19M=2.8'X+14.-2./X 
IF(X.LT.,5)  M=11.8+ALOG10(X) 
IF(M.LT.M4)  M-M4 
BJ(M)«0. 

JM-M-1 

BJ(JM)x1. 

DO  16  J-4.M 
J2-JM 
JM-JM-1 
J1-JM-1 

BJ(JM)-J  1  /X'BJ(J2)-BJ(JM'f2) 

16  CONTINUE 
SsO. 

IF(M.LE.4)  GOTO 24 
DO  17  J=4.M,2 
S-S4,BJ(J) 

17  CONTINUE 

24  S»BJ(2)+2.*S 

18  ARG-Z2(L)'CC 


UA=AT(L)/S"CMPU(COS(ARG).SIN(ARG)) 
UB-XT(L)-UA 
DO  25  MaM3.M4 
FA(M)-BJ(M)'UA«FA(M) 

FB(M)-BJ(M)'UB4^FB(M) 

25  CONTINUE 
15  CONTINUE 

IF(M1.NE.0)  GO  TO  26 
FA(1)-  FA(3) 

FB(1)--  FB(3) 

26  UA-U1 

DO  27  M-M5.M6 

M7-M-1 

M8-M+1 

F2A>UA’<FA(M8)«FA(M7)) 

F2B-UA‘(FB(M8)+FB(M7)) 

UB-U*UA 

F1A-UB’(FA(M8)-FA(M7)) 

F1B>UB*(FB(M8)-FB(M7)) 

U4<D4’UA 

U2-D3*F1A+U4‘FA(M) 

U3-D3*F1B+U4-FB(M) 

U4=DR*F2A 

U5.0R*F2B 

K1-K2-1 

K4-KUN 

K5-K2+N 

R(K2*MT)=-D5*(F2A+D2'F2B) 

R(K5+MT)-Dr  (F 1  A+D2*F  1 B) 

IF(IP.EQ.1)G0T0  21 
R(K1)-R(K1)4^U2-U3 
R(K4)=R(K4)+U4-U5 
IF(IP.EQ.MP)  GO  TO  22 

21  R(K2)=U2+U3 
R(KS)>U4+U5 

22  K2-K2+N2 
UA-UB 

27  CONTINUE 
12  CONTINUE 

RETURN 

END 

SUBROUTINE 

PLANEW(NWIRES,NW1  .NW2.NP1  ,NP2.XH.YHZH, 
•THR.PHR.R) 

C 

C  PLANE  WAVE  EXCITATION  VECTOR  ELEMENTS  FOR 
WIRE  AND 

C  ATTACHMENT  SEGMENT,  INCIDENCE  DIRECTION  IS 
(THR.PHR). 

C  -  ONLY  ONE  INCIDENCE  ANGLE  PER  CALL 

C  -  WIRE  DOES  NOT  HAVE  TO  LIE  IN  THE  Z-0  PLANE 

C 

COMPLEX 

U0.AA.BB,C.R(400).CEXP,EXP,FI1,FI2.SI,DI,CMPLX 

DIMENSION 

ZH(400).NS(2),NW1(2).NW2(2),XH(400),YH(400) 
MP2-NP2-1 
MT2-NP2-2 
MT22=2*MT2 
DOSL-1.NWIRES 
5  NS(L)-NW2(L)-NW1(L)f1 
U0-(0..0.) 

CC-COS(THR) 

SS-SIN(THR) 

CP-COS(PHR) 
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SP>SIN|PHR) 

UP.SS*CP 

VP.SS’SP 

WP«CC 

DO  12  IP=1,MP2 

ll«IP+NP1 

I-II+1 

ZS..5-(ZH(I)+ZH(11)) 

XS«.S*(XH(I)+XH(M)) 

YS«.S*(YH(I)+YH(II)) 

OX>XH(l)-XH(ll) 

DY«YH(I)-YH(II) 

D1.SQRT(DX"2+DY"2) 

SU=DY/D1 

CU=DX/D1 

C  FOR  WIRES  IN  PERPENDICULAR  TO  Z  SIN(V)=1  AND 
COS(V)»0 
SVsI.O 
CV=0.0 
C 

C  WIRE  SEGMENT  CALCULATIONS 
C 

A.UP*CU+VP’SU 

B»UP*XS+VP'YS+WP*ZS 

C=CMPLX(0.A) 

EXP=CEXP(CMPLX(0.,B)) 

AA=CC*(CU'CP+SU-SP) 

BB=SU*CP-SP-CU 

PSI=DrA2. 

IF(PSI  NE.O.)  GO  TO  60 
SINC=1, 

GO  TO  61 

60  SINC*SIN(PSI)'PSI 

61  COSP=COS(PSI) 

FI1=SINC’DrEXP/2. 

FI2=(0..0.| 

IF(ABS(A).LT.1.E-4)  GO  TO  62 
CSP=COSP-SINC 
IF(ABS(CSP).LT.1.E-4)  GO  TO  62 
FI2=EXP'C-CSP 

62  CONTINUE 
SI=FI1+FI2 
DI=FI1-FI2 

C 

C  R-WIRE-THETA 
C 

IF(IP.EQ.MP2)  GOTO  10 
R(1P)=AA‘SI 
R(IP+MT2)=BB'SI 
10  CONTINUE 
C 

C  R-WIRE-PHI 
C 

14IF(IP.EQ.1)G0T0  12 
R(IP-1)=R(IP-1)+AA'DI 
R(IP-1+MT2)=R(IP-1+MT2)+BB‘DI 
12  CONTINUE 
RETURN 
END 

Raflactor 

PROGRAM  mrc.f  (Marin*  Radio  Ant«nna  Program) 
c  Modifiod  from  ORDONEZ.f  -  hyperboloid  baconws  a  disc 


C  Paraboloia  WITH  A  PARASITIC  RING  DIPOLE  FEED 
Fad  dpole  a  at  z=0. 

C  The  program  usas  a  parasitic  dipola  and  two  h-plan* 

C  parasitic  dipolas  and  the  subraflector  This  system  also 

contains 

C  a  cavity  nng. 

C 

C  ICALCaO  CALC  CURRENTS  AND  FIELD 
C  >1  CALC  FIELDS  ONLY  (READ  CURRENTS) 

C  IMP-0  PEREECTLY  CONDUCTING  SURFACES 
C  -1  SOME  NONZERO  SURFACE  RESISTANCE 
C  ICWRT-0  WRITE  CURRENT  COEFFICIENCTS  TO 
FILE  peurrant 

C  IRES-0  READ  RESISTIVE  CORRECTIONS 
C 

COMPLEX 

EP.ET^(500000),RS(1000).RW(400).E(700).C(700).U.UC 
COMPLEX  ZL0(400)ZL(1500).EC.EXZSS(SOO0O| 
COMPLEX  EXP1.EXP22SW(S0OOO)ZWW(10OO0) 
COMPLEX 

CEXP.CONJG.CMPLX.ET1  ,ET2.EP1  .EP2.ETW.EPW 
DIMENSION 

RH(400)ZH(400).XT(4)AT(4).IPS(700).NW1(4).NW2(4) 

DIMENSION 

A(400).X(400).EXP(800)ANG(800).ECP(800) 

DIMENSION  ECV(800).EXV(800).PHC(800).PHX(800) 
DIMENSION 

XX(41).YY(41).XH(400).YH(400).PHA(400) 

DATA 

PI.START.STOP.DT,3.141S9.0..180..1./,ICWRT  0  .IRES  1 
data  ICALC-'O'.IPRINT  OMMP/O'.ITEST.T. 

DATA  nVZ',xt(  1).at(1)'.5773S03.1 
c  data  nt4,  .1(1(1  ).irt(2).at(1),at(2)' 33998104,861 13631 15, 
c  •  .6521451548.3478548451' 

OPEN(1,FILE-'outgaus'.status-'old') 

OPEN(30.FILE»'rh.m) 

OPEN(31.FILE-zh.m) 

READ(1.*)  NPHI 
D0  3K«1,NPHI 
READ(1.')X(K).A(K) 

3  CONTINUE 
RAD-PIT80. 

BK-2.*PI 

U=(0.1.) 

U0=(0..0.) 

UC— U/4./PI 
NT2-NT/2 
DO  1  K=1,NT2 
K1-NT-K+1 
AT(K1).AT(K) 

XT(K1)=XT(K) 

1  XT(K)--XT(K) 

C 

C  NP10  -  NUMBER  OF  PARABOLOID  POINTS 
C  NP1 1  •  NUMBER  OF  DISC  POINTS 
C  NP12  -  NUMBER  OF  RING  POINTS 
C 

MODES-1 

MHI-MODESfl 

NBLOCK-2*MODES+1 

NWIRES-1 

WRITE(6.*)  ’ENTER  START.STOP.DT’ 

READ(5.’)  START.STOP.DT 

NP10-61 

NP11-31 
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NP12-3 

WRITE(6.*)  ENTER  PARA  (NP10),  DISC  (NP1 1)  AND 
RING  (NP12)  PTS' 

READ(S,*)  NP10,NP11.NP12 
NP1-NP10+NP11+NP12 
C 

C  DM«MAIN  REFLECTOR  DIAMETER 
C  DS-DISC  (SUBREFLECTOR)  DIAMETER 
C  ZCsDISC  LOCATION  ALONG  Z  AXIS 
C  FOD«f/DM  RATIO 
C  RC-RADIUS  OF  THE  RING 
C  SW«RING  (STRIP)  WIDTH 
C 

DM>5.0 

R0>DM/2. 

DS=1,1 

ZC«.25 

FOD=.3 

RC=0.5 

SW-.05 

WRITE(6.-)  'ENTER  DM.  DS,  ZC  (>0).  FOD,  RC.  ZR 
AND  SW 

READ(5.*)  DM,DS2C,FOD.RC.ZR,SW 
FM=FOD‘DM 

WRITE(6.’)  'ENTER  ZSHIFT  FOR  PARABOLOID' 
READ(5,*)  ZSHIFT 
C 

C  GENERATE  PARABOLOID  CONTOUR 
C 

PHIV=2.*ATAN(1  ..^./FOD) 

FNM=FLOAT(NP10-1) 

DO  50  l=1.NP10 

TH=FLOAT(l-1)‘PHIV/FNM 

RM.2.‘FM/(1.+COSfrH)) 

XH(l)=0. 

YH(()=0. 

ZH(l)=-RM'COS(TH)*BK+ZSHIFT*BK 
RH(I)=RM*SIN(TH)"BK 
PHA(ll)=0 
50  CONTINUE 

DISC  LOCATED  A  Z=ZC 

DSK=DS/2./FLOAT(NP1 1  - 1 ) 

D0  5I=1.NP11 
llsl+NP10 
XH(ll)=0. 
yH(ll)=0. 

ZH(II)=ZC'BK 

RH(  ll)=FLOAT(  I- 1  )*DSK’BK 
PHA(ll)=0. 

5  CONTINUE 
C 

C  NEXT  THE  SIDEWALLS 
C 

DSW=SW/FLOAT(NP1 2- 1 ) 

ZRBK=ZR*BK 
IF(NP12.NE.O)  THEN 
DO?  I=1,NP12 
II.NP11+NP10+I 
XH(ll).0. 

YH(ll)=0. 

RH(II)«RC*BK 

ZH(ll)«FLOAT(l-1)*DSWBK+ZRBK 
7  PHA(ll).0. 


ENDIF 

NP1-NP10+NP11+NP12 

C 

C  DIPOLE  ANTENNA  POINTS  IN  WAVELENGTHS:  AW. 
RADIUS:  NP2.  NUMBER 

C  OF  POINTS  (.NP2UNP22:  NP21  MUST  BE  ODD); 
DLGTH.  LENGTH; 

C  PLGTH  IS  PARASITIC  ELEMENT  LENGTH 
C 

AW..001 
OLGTH..45 
PLGTH.  44 

WRITE(6,')  ENTER  DLN.  PlN.  DDIP,  ROT  AND  AW 
WRITEIS,*)  '(DDIP<0  IS  TOWARD  PARABOLOID)' 
READ(S.‘)  DLGTH.PLGTH.DDIP,R0T.AW 
CENT=DLGTH,'2. 

AK.AW*BK 

NP21.il 

NP22>K3 

WRITE(6,')  ENTER  NP21  AND  NP22  (ONLY  1  DIP  IF 
NP22.0) 

READ(5.*)  NP21.NP22 
NW1(1)=NPU1 
NW2(1)=NP1+NP21 
NP2.NP2UNP22 
DEL.DLGTH/FLOAT(NP21  ■  1 ) 

DO  10  I.1.NP21 

II.NPUI 

YH(ll).0. 

ZH(ll)=0 

XH(ll)=BK‘(FLOAT(l- 1  )*DELCENT) 
RH(II).SQRT(XH(II)”2+YH(II)**2) 
PHA(II)=ATAN2(YH(II).XH(II)+1  E-6) 

10  CONTINUE 
IF(NP22.NE.O)  THEN 
NWIRES.2 
NW1(2).NP1+NP21+1 
NW2(2).NP1  +NP2 1 +NP22 
DEL=PLGTHELOAT(NP22- 1 ) 

CENT.PLGTH2. 

DO  11  I.1.NP22 
II.NP1+NP21+I 
ZH(II)=DDIP'BK 

XH(ll).BK‘(FLOAT(l-1)-DEL-CENT)*COS(ROrRAD) 

YH(ll)=BK’(FLOAT(l-1)’DEL-CENT)*SIN(ROT’RAD) 

RH(ll).SQRT(XH(ll)'*2+YH(li)"2) 

PHA(II)-ATAN2(YH(II),XH(II)+1.E-6) 

11  CONTINUE 
ENDIF 

C  NP23.NUMBER  OF  H-PLANE  DIPOLE  POINTS 
C  HLGTH.H-PLANE  DIPOLE  LENGTH 
C  DHP.H-PLANE  DIPOLE  DISTANCE  FROM  FEED 
C  ZHP.H-PLANE  DIPOLE  Z  COORDINATE 
C  NP24.NUMBER  OF  POINTS  FOR  THE  OTHER 
H-PLANE  DIPOLE 
NP23-7 
NP24.7 
HLGTH..2 
DHP..2 
ZHP..001 

WRITEie.*)  ENTER  NUMBER  OF  POINTS  PER 
H-PLANE  DIPOLE' 

WRrrE(6.*)  '(0  FOR  NO  H-PLANE  PARASITIC 
DIPOLES)' 

READ(5.-)  NP23 
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IF(NP23.NE.O)  THEN 

WRITE(6.‘)  H-PLANE  DIPOLE  LENGTH.  DISTANCE 
AND  LOCATION' 

REAO(S,*)  HLGTH  DHP^HP 

NWIRES>NWIRES«2 

NP24.NP23 

NP2.NP21+NP22+NP23+NP24 
NW1(3)»NPUNP21+NP22+1 
NW2(3)-NP1  +NP2 1  ♦NP22+NP23 
NW1  (4).NP1  +NP21+NP22+NP23+1 
NW2(4).NP1+NP2>*NP22+NP23+NP24 
DEL»HLGTH/FLOAT(NP23- 1 ) 

CENT=HLGTH,2. 

DO  12I-1.NP23 
IUNP1+NP21+NP22+I 
ZH(l1)sZHP*BK 
YH(I1)«DHP'BK 

XH(  1 1  )«BK*(FLOAT(l- 1  )*DEL-CENT) 
RH(l1)=SORT(XH(l1)“2+YH(l1)"2) 

PHA(  1 1  )»ATAN2(  YH(I1 )  ,XH(  II  )♦  1  .E-6) 

I2=IUNP23 

ZH(I2)«ZHP*BK 

YH(I2)=-DHP*BK 

XH(l2)=BK'(FLOAT(l-1)‘DEL-CENT) 

RH(I2)=SQRT(XH(I2)"2+YH(I2)"2) 

PHA(l2)=ATAN2(YH(l2).XH(l2)+1.E-6) 

12  CONTINUE 
ENDIF 
CLOSE!  1) 

NP=NP1+NP2 

WRITE(6.*)  ENTER  FREQUENCY  FACTOR' 

READ(5.’)  FFAC 
OPEN(8.fil«='outmrc') 

WRITE(8,8000) 

DM.DS.FOD.PHIVWAD.NP10.NP11.RC.SW, 

NP12ZSHIFT,OLGTH,PLGTH.ODIP.NP2,AW,MODES.NT. 

NPHI.FFAC 

IF(NP23.NE.0)  WRITE(8.8001)  HLGTH.DHPZHP,NP23 

8000  FORMAT('/,SX.'’"  MARINE  RADIO  REFLECTOR 
ANTENNA  5X. 

"REFLECTOR  PARAMETERS  (LENGTHS  IN 
WAVELENGTHS):', /.5X. 

"MAIN  REFL  DIA='.F8.3.'.5X.'DISC 
DIA='.F8.3./,5X.'F/0='.F8.4. 

•  /.SX.'PHIV  (DEG)='.F8.2./.5X.'MAIN  REFL  POINTS 
(NP11)=',I4./, 

•  SX.'SUB  REFL  POINTS  (NP12)=',I4,  ,5X./.5X, 

■  'RING  RADIUS='.F8.3.'.5X, 'SIDEWALL 

LENGTH=',F8.3./.5X. 

•  'NUMBER  OF  RING  POINTS  (NP12)='.I4.//.5X. 

•  PARABOLOID  SHIFT  DISTANCE. 

ZSHIFT=',F8.3./.5X, 

•  DIPOLE  PARAMETERS:',/.5X.  FED  DIPOLE 
LENGTH-', F8.3./.5X, 

•  PARASITIC  DIPOLE 
LENGTH-'.F8.3./.SX,'SPACING-',F8.3./.5X. 

•  'NUMBER  OF  DIPOLE  POINTS  -',I4,/,5X.  RADIUS  OF 
DIPOLE-', 

•  F8.4.//.5X.'NUMBER  OF  AZIMUTHAL 
MODES.'.l3./,SX.'NT='.l3. 

•  /.5X.'NPHI.'.I3,/,5X.'FREQ.  SCALE  FACTOR-’,F8.4) 

8001  FORMAT(/.5X.'H-PLANE  DIPOLE 
LENGTH='.F8.3,/,5X. 


•  'H-PLANE  DIPOLE  DISTANCE  FROM 
AXIS- .F8.3./.5X. 

•  Z  POSITION  OF  H-FIELD 
DIPOLE-'.F8.3.'.SX .'NUMBER  OF  H-FIELD 

•  DIPOLE  POINTS-  ,I3) 

IF(ISEQ.EQ.O)  WRrrE(8.1300) 

1300 

FORMAT!/ .5X,'INDEX',8X,7(I)'.10X.'RHO(I).8x.ZSURF) 

c 

C  FILE  dMat  CONTAINS  RESISTIVE  CORRECTIONS  IF 

IRESNEO 

C 

IF((IMP.EQ.O).OR.(IRES.NE.O))  GO  TO  70 
OPEN(9.FILE-'iMdat') 

READ(9.‘)  ID 
WRITEIS.-)  ID-'.ID 
DO  75  1-1, ID 

READ(9,*)  II.RI.YY(I),DUM 
IF(YY(I).LT.0.)  YY(I).0. 

75  CONTINUE 

XXI-I./FLOATIID-I) 

DO  15  11-1. ID 
XX(ll)-XXrFLOAT(ll-1) 

IS  CONTINUE 

WRITE(6,-)  SURFACE  IMPEDANCE  VALUES  READ 
70  CONTINUE 
DO  52  1-1  .NP1 

IF(ABS(ZH(I)).LT..001)  ZH(l)=0. 
IF(ABS(RH(I)),LT..001)  RH(l)ri) 

ZLO(I)-(0„0.) 

XH(I)-XH(I)*FFAC 
YH(I)-YH(I)*FFAC 
ZH(I)-ZH(I)*FFAC 
RH(I)=RH(I)’FFAC 
ZHB-ZH(I)BK 
RHB»RH(I)/BK 
WRrrE(30.5019)  RHB 
WRITE(31.5019)ZHB 
c  if(rt*.9».10  )  Il0(i)»(.2.0.) 

C 

C  NONZERO  SURFACE  IMPEDANCES  IF  IMP  >  0 
C 

IF((IMP.EQ.0)  OR.(IRES.NE.O))  GO  TO  51 
IF(I.GT.NPIO)  GO  TO  51 
C 

C  USE  THE  EXACT  RESISTIVE  CORRECTIONS  - 
INTERPOLATE  BETWEEN 
C  DATA  POINTS. 

C 

UU-RHB/RO 

CALL  INTERP(UU.VV,XX.YY,ID) 

ZL0(l)=VV 

51  CONTINUE 

WRITE(8,e004)  IZHB.RHBZLO(I) 

52  CONTINUE 

8004  F0RMAT{6X.I4.4X,F8.3,8X.F8.3,6X,2F8.4) 
WRrTE(8.1310) 

1310  F0RMAT(/,5X,’DIP0LE  COORDINATES:',//, 

•  5X,'INDEX',8X,Z(I)',10X,'RHO(I)',10X,'PHI  PLN  ) 

DO  53  l-NP1+1,NP 
IF(ABS(ZH(I)).LT„001)  ZH(l)-0. 
IF(ABS(RH(I)).LT.,001)  RH(I)-0, 

XH(I)=XH(I)’FFAC 

YH(I)-YH(I)*FFAC 

ZH(I)-ZH(I)TFAC 
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RH(I)«RH(I)TFAC 
PHB»PHA(I)/RAD 
ZHB«ZH(I)/BK 
RHB«RH(I)/BK 
wRrrE(30.soi9)  rhb 
WRrrE(31.S019)  ZHB 
53  WRrrE(8.B00S)  UHB.RHB.PHB 
800SFORMAT(6X,I4.4X,P8.3.8X.F8.3.9X.F8.2) 

799  CONTINUE 
IF(ITEST.EQ.O)  GO  TO  9998 
wnt«(6,‘)  gaomatry  dafinad' 

C 

C  DEFINE  DIMENSIONS  OF  THE  IMPEDANCE  MATRIX 

BLOCKS 

C 

MT1-NP1-2 

MP1-NP1-1 

N-MP1+MT1 

MT2.NP2-2 

NSURF>NBLOCK*N 

NSW=NSURF+MT2 

NROW«NSW 

WRITE(6.‘) 

NPI  ,NP2,NSURF,NROW='.NP1  .NP2.NSURF.NROW 
WRITE(6,*)  'MT2  SHOULD  BE  ODD  -  MT2='.MT2 
IF(ITEST.EQ.O)  GO  TO  9998 
IF(ICALC.EQ.O)  THEN 
C 

C  COMPUTE  IMPEDANCE  MATRIX  ELEMENTS 
C 

CALL 

ZMATWW(NWIRES.NW1,NW2,XH,YH,RHZH,NT.XT.AT.A 

K.ZWW) 

CALL  ZASMB0(NP1  ,NP2,MODES7.ZWW) 
wrila(6.*)  'wira  impadanca  computad' 

IF(IMP.NE.O)  THEN 

CALL  ZLOAD(NP1  ,RH2H2L0ZL) 

WRrTE(6.')  RESISTIVE  SURFACE  IMPEDANCE 
MATRIX  COMPUTED' 

ENDIF 

DO  400  M=1,MHI 

NMsM-1 

CALL 

ZMATSS(NM.NM,NP10,NP11,NP12.NPHI.NT,RHZH.X.A,X 

T.AT,ZSS) 

IF(IMP.NE.O)  CALL  ZTOTIMTI  .MP1 2L2SS) 

CALL 

ZMATSW(NWIRES.NW1,NW2,NP10.NP11.NP12.XH.YH.R 

H2H, 

•  NT,XT,AT,NPHI.X,A,NM.AK2SW) 

CALL  ZASMBN(NP1  ,NP2,MODES.NM2ZSS2SW) 

IF(NM.EQ.O)  GO  TO  400 

NMN>-NM 

CALL 

ZMATSW(NWIRES,NW1,NW2.NP10,NP11.NP12,XH.YH.R 

H2H, 

•  NT.XT.AT.NPHI,X.A,NMN.AK2SW) 

CALL  ZASMBNINPI  ,NP2.MODES.NMN22SS2SW) 
400  CONTINUE 
wnta(6.‘)  'z  fillad' 

CALL  DECOMP(NROW.IPS2) 
wrta(6,*)  'z  dacompoaad’ 

C 

C  EXCITATION  ELEMENTS  (NONZERO  FOR  THE 
DIPOLE  ONLY) 


C 

DO40  U1.NROW 
40  B(l)>(0..0.) 

IO>rfJSURF+(NP21t)2 

B(l0)>(1..0.) 

CALL  SOLVE(NROW.IPS2.B.C) 
writa(6,‘)  'currant  coafRcianls  computad' 
c  do  665  iiBl.ivow 
c665  writa(6,*)  ii.clii) 

IF((ICVmT.EQ.0).AND.(ICALC.EQ.O))THEN 

C 

C  WRfTH  CURRENTS  ON  DISC  FOR  PATTERN 
INTEGRATION 
c  (program  cavarchint.r) 

C 

OPEN(3.NaB'pcurrant') 

WRITE(3.*) 

OM^W2C,RC,DS.NP1,NP2.NBLOCK,NWIRES.FFAC 

VYRrTE(3.') 

NP.MT1  .MP1  .N.MT2.NSURF.NROW.MHI.MODES 
VYRITE(3.')  (XH(I).I«1,NP) 

VYRrrE(3,*)  (YH(I).I.1,NP) 

VVRITE(3,*)  (ZH(I).I.1.NP) 

WRrrE(3,*)  (RH(I).I=1.NP) 

WRITE(3.*)  (PHA(I),I-1,NP) 

VVRITE(3.*)  (NW1(I),I=1.NWIRES) 

WRITE(3.')  (NW2(I),I»1,NWIRES) 

WRITE(3.*)  (C(l).l=1.NROW) 

ENDIF 

ENDIF 

C 

C  IF  ICALC.NE.O  THEN  READ  CURRENT  COEFFICIENTS 
FROM  DISK  FILE 
C 

IF(ICALC.NE.O)  THEN 
OPEN(3.fila«'pcurTant',STATUS«'OLD') 

READ(3.*) 

DMX.AWX  2CX.RCX.DSX.NP  1 X  .NP2X  .NBLX.NWX 
READ(3,*) 

NPX.MT1X.MP1X.NX,MT2X,NSX.NRX,MHIX.MODESX 
READ(3.’)  (xh(i),i.1,rp) 

READ(3,-)  (yh(j).i«1.np) 

READ(3.*)  (zh(i).i=1.rp) 

READ(3.’)  (rh(i).bl.np) 

READ(3.*)  (pha(i).i»1.np) 

READ(3.’)  (nwl(i).bl.nwiras) 

READ(3.*)  (nw2(i),i=1.nwiras) 

READ(3.")  (clil.ttt.nrow) 

CLOSE(3) 

writa(6.*)  'data  raad  from  pcurrant' 
do  909  bl.np 

909  writa(6.')  'i.rh.zhs',i.rt)(i)/bk.zh(i)/bk 
ENDIF 

IT«INT((STOP-START)/DT)+1 

C 

C  RECEIVER  PHI  CUTS;  DO  0  AND  90  (E-  AND  H- 
PLANES) 

C 

DO  501  IP«1,2 

ECX>1.a-10 

PHRO-0. 

IF(IP.EQ.2)  PHR0«RAD*90. 

DO  500  l«1.rT 

THETA.FLOAT(l-1)*DT+START 

THR»THETA*RAD 
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PHR-PHRO 

IFfTHETA  LE  ISO.)  GO  TO  99 
THR.(360 -THETA)*RAD 
PHR.PHRO+PI 
99  CONTINUE 

c  wnt*(6,*)  Ihata.phK'.thata.phr/rad 
ET1-U0 
ET2-U0 
EPUUO 
EP2-U0 
ETW-UO 
EPW=U0 

DIPOLE  FIELD 

CALL 

PLANEW(NWIRES,NW1  ,NW2,NP1  ,NP2,XH,YH^H.THR.P 
HR.RW) 

DO210L»1,MT2 

ETW.ETW+RW(L)*C(L+NSURF) 

210  EPW=EPW+RW(L+MT2)*C(L+NSURF) 

C 

C  REFLECTOR  AND  CAVITY  FIELD  CONTRIBUTIONS 
C 

00  300M=1.MHI 
NM=M-1 

EXP1=CEXP(CMPU(0.,FLOAT(NM)‘PHR)) 

EXP2»CONJG(EXP1) 

CALL 

PLANES(NM,NP10.NP11,NP12,NT.RH2H,XT,AT.THR.RS) 
NTOP1.MODES-NM 
NTOP2=NBLOCK-(NTOP  1  ♦  1 ) 

NS2=NTOPrN 
NS1=NTOP2'N 
DO250L«1,MT1 
ET1  aETI  +RS(L)*C(L+NS1  )-EXP  1 
EP1  «EP1  ♦RS(L+N)’C(L+NS1  )’EXP  1 
IF(NM.EQ.O)  GO  TO  2S0 
ET1=ET1+RS(L)’C(L+NS2)’EXP2 
EP1=EP1-RS{L+N)'C(L+NS2)'EXP2 
250  CONTINUE 
00  260L=1.MP1 

ET2»ET2+RS(L+MT1)‘C(L+NSUMT1)*EXP1 
EP2=EP2+RS(L+MT1+N)*C(L+NS1+MT1)*EXP1 
IF(NM.EQ.O)  GO  TO  260 
ET2»ET2-RS(L+MT1  )*C(L+NS2+MT1  )’EXP2 
EP2=EP2+RS(L+MT1  +N)’C(L+NS2+MT1  )‘EXP2 
260  CONTINUE 
300  CONTINUE 

ET-UC*(ET1  +ET2+ETW) 

EP=UC'{EPUEP2+EPW) 

EC=ET 

EX=EP 

ECV(I)=CABS(EC) 

EXV(l)sCABS(EX) 

ECR.REAL(EC) 

ECI-AIMAG(EC) 

EXR»REAL(EX) 

EXI-AIMAG(EX) 

PHC{I)»ATAN2(ECI,ECR+1  .•-10)/RAD 
PHX{I)«ATAN2(EXI,EXR+1.«-10)/RAD 
ANG(I)»THETA 

ECX»AMAX1(ECX,ECV(I),EXV(I)) 

500  CONTINUE 
WRrrE(8.103)  PHR/RAD,ECX 


103  FORMAT(/,5X,  ’PHI  OF  RECEIVER 
(OEG)-'.F8.2./.5X. 

*  MAXIMUM  FIELD  VALUE  (V/M).  £15.5) 

D0  600  l«1.IT 
ECP(I).(ECV(I)/ECX)"2 
EXP(I)»(EXV(I)/ECX)"2 
ECP(I)>AMAX1(ECP(I)..00001) 

EXP(I)>AMAX  1  (EXP(I) .  .00001 ) 

ECP(I)>10.‘ALOG10(ECP(I)) 

EXP(I)«10.‘ALOG10(EXP(I)) 

600  CONTINUE 
IF(IPRINT.NE.O)  GO  TO  310 
WRITE(8.5015) 

5015 

FORMAT(///.7X  .'ANGLE'.  1 5X  ,'CO-POL'.25X  ,X-POL'./.7X , 
•'(DEG)',4X,2('(VOLTS)',4X,'(DEG)'.3X,'(DB-REL)',4X)) 
D0  9000L«1.rT 
VYRrrE(8.5016) 

ANG(L).ECV(L),PHC(L),ECP(L).EXV(L).PHX(L) 

•.EXP(L) 

5016FORMAT(SX.F6.1.3X.2(F8.4,3X.F7.1.3X.F7.2.3X)) 
9000  CONTINUE 
310  CONTINUE 
C 

C  WRITE  MATLAB  FILES  (  .m  FILES) 

C 

IF(IP.EQ.1)THEN 

OPEN|2.lilc-'ang.m') 

OPEN(3,lil*B'q>ol*.m') 

OPEN(4.fil*«'xpol*.m') 

DO  9097  l«1.IT 
WRITE(2.S019)  ANG(I) 

WRITE(3.5019)  ECP(I) 

9097  WRITE(4,5019)  EXP(I) 

CLOSE(2) 

CLOSE(3) 

CLOSE(4) 

ENDIF 

IF(IP.EQ.2)  THEN 
OPEN(3.fil««'cpolh.m') 

OPEN(4,fil»«'xpolh.m') 

DO  9098  1=1. IT 
WRITE(3.5019)  ECP(I) 

9098  WRITE(4,5019)  EXP(I) 

CLOSE(3) 

CLOSE(4) 

ENDIF 

501  CONTINUE 
5019  FORMAT(F8.3) 

9998  STOP 
END 

SUBROUTINE  SOLVE(N,lPS,UL,B,X) 

COMPLEX  UL(500000).B(700),X(700).SUM 
DIMENSION  IPS(700) 

NP1=N+1 

IP.IPS(1) 

X(1)=B(IP) 

DO  2  l=2.N 
IP.IPS(I) 

IPB«IP 

IM1>I-1 

SUM<(0..0.) 

DO  1  J>1.IM1 
SUM»SUM+UL(IP)*X(J) 

1  IP-IP+N 
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2  X(I)>B(IPB)-SUM 
K2-N’(N-1) 

IP>IPS(N)-t'K2 

X(N)«X(N)/UL(IP) 

OO  4  IBACK«2.N 

l>NP1-IBACK 

K2>K2-N 

IPI>IPS(l)-fK2 

IP1-I+1 

SUM>(0.,0.) 

IP>IPI 

D0  3J-IP1.N 
IPsIP-fN 

3  SUM>SUM+UL(IP)'X(J) 

4  X(I)>(X(I)-SUM)/UL(IPI) 

RETURN 

END 

SUBROUTINE  DECOMP(N.IPS,UL) 

COMPLEX  UL(500000), PIVOT, EM 
DIMENSION  SCL(700).IPS(700) 

D0  5U1.N 

IPS(I)=I 

RNsO. 

Jlsl 

D0  2J«1,N 

ULM>ABS(REAL(UL(J1 )))+ ABS(AIMAG(UL(J1 ))) 

J1>JUN 

IF(RN-ULM)  1,2.2 

1  RN>ULM 

2  CONTINUE 
SCL(I)=1./RN 

5  CONTINUE 
NM1«N-1 
K2x0 

D0  17K=1,NM1 
BIGxO. 

DO  11  IxK.N 
IP-IPS(I) 

IPK»IP+K2 

SIZE=(ABS(REAL(UL(IPK)))+ABS(AIMAG(UL(IPK))))’SCL( 

IP) 

IF(SIZE-BIG)  11,11,10 

10  BIGxSIZE 
IPV=I 

1 1  CONTINUE 
IF(IPV-K)  14,15,14 

14  JxlPS(K) 

IPS(K)=IPS(IPV) 

IPS{IPV)«J 

15  KPP«IPS(K)+.K2 
PIVOT«UL(KPP) 

KP1-K+1 

DO  16  l»KP1,N 
KP-KPP 
IP>IPS(I)4'K2 
EM»-UL(IP)/PIVOT 
18UL(IP)«-EM 
DO  16  J«KP1,N 
IP-IP+N 
KP«KP-fN 

UL(IP)-UL(IP)+EM*UL(KP) 

16  CONTINUE 
K2xK2-fN 

17  CONTINUE 


RETURN 

END 

FUNCTION  BLOG(X) 

IF(X,GT..1)GOTO  1 
X2«X'X 

BLOG>((.075*X2-.1666667)-X2'f1  .)'X 
RETURN 

1  BL0G>AL0G(X4GQRT(1>X’X)) 

RETURN 

END 

SUBROUTINE  2LOAD(NP,RH2HZOZ) 

C 

C  COMPUTES  IMPEDANCE  MATRIX  ELEMENTS  FOR 
LOADED  BODIES  OF  REV 
C  Z0(l)  IS  THE  SURF  IMPEDANCE  OF  THE  ITH 
FOMENT  (NP-1  SEGMENTS) 

C  Z(.)  ARE  THE  IMPEDANCE  MATRIX  TERMS 
{TRIDIAGONAL  FOR  T-T 

C  SUBMATRIX;  DIAGONAL  FOR  P-P  SUBMATRIX) 
STORED  IN  COL  VECTOR 
C 

COMPLEX 

C1,C2iO{400)2(1500)J(1,X2,X3,Y1.y2,Y3,FN(400) 
COMPLEX  U1,U2,U3,XI,YI 

DIMENSION  RH(400)2H(400),RS(400).D(400).SV(400) 

Pl>3.14159 

MT.NP-2 

MP«NP-1 

N>MT+MP 

DO  10  IP«2,NP 

llxlP-1 

Dfl»RH(IP)-RH(ll) 

DZ«ZH(IP)-ZH(II) 

D(II)»SQRT(DR*DR+DZ*DZ) 

SV(II)«DR/D(II) 

RS(ll)«.S‘(RH(IP)-i'RH(ll)) 

DS.D(ll)’SV(ll)/2. 

Q1>RS(II)+DS 

Q2>RS(II)-DS 

FN(II)«1 

IF((ABS(Q2).GT.1.E-6).AND.(ABS(Q1).GT.1.E-6)) 

*  FN(ll)>ALOG(Q1/Q2) 

10  CONTINUE 
L0-MT*3-2 
DO20l>1.MP 
CIxPI'ZOd) 

IF(I.EQ.MP)  GO  TO  80 
Klx2 

IF(I.EQ.1)KI.1 
IF(I.EQ.MT)  KI.3 
lbl-^1 

C2.PrZ0(ll) 

A»SV(I) 

IF(ABS(A).LT.1  .E-6)  GO  TO  41 

X1.C1*FN(l)/2./A 

X2«C1*2,/A*(1.-RS(I)*FN(I)/D(I)/A) 

X3«-X2’RS(I)ID(I)/A 
GO  TO  42 

41  CONTINUE 
X1«C1/2./RS(l)*0(l) 

X2>(0..0.) 

X3«CrD(l)/6./RS(l) 

42  CONTINUE 
A>SV(II) 

IF(ABS(A).LT.1.E-6)  GO  TO  45 
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Y1«C2‘FN(ll)/2,/A 
Y2.C2‘2./A*(1.-RS(II)*FN(IIMD{II)/A) 
Y3«-Y2*RS(II)/D(1I)/A 
GO  TO  40 
45  CONTINUE 
Y1«C2«./RS(II)’0(II) 

Y2>(0..0.) 

Y3«C2*D(ll)/6./RS(ll) 

40  CONTINUE 
C 

C  DEFINE  TRIDIAGONAL  ELEMENTS  FOR  T-T 
SUBMATRIX  (STORED  IN  COLS) 

C  (U1-  DIAG;  U2-  LOWER;  U3-  UPPER) 

C 

XI>XUX2+X3 

YI.Y1-Y2+Y3 

IF(KI.EQ.1)XI«C1/SV(I) 

C  IF(KI.EQ.3)  YI=C2/SV(II) 

UlxXI+YI 

U2»X1-X3 

U3«Y1-Y3 

L«2+(l-2)-3 

IF(KI.EQ.I)  LxO 

LUL+1 

L2-L+2 

L3xL+3 

go  to  (50,50.70), ki 
SO  2(L1)=U1 
Z(L2)=U2 
GO  TO  80 
60Z(L1)xU3 
Z(L2)=U1 
Z(L3)xU2 
GO  TO  80 
70Z(L1)«U3 
Z(L2)«U1 

80  Z(L0*I)=2  •CrO(l)/RS(l) 

20  CONTINUE 
RETURN 
END 

SUBROUTINE  ZTOT(MT,MP.ZL2) 

C 

C  ADDS  THE  SURF  IMPEDANCE  TERMS  TO  THE 
TRIDIAGONAL  ELEMENTS  OF 
C  THE  BOR  IMPEDANCE  MATRIX  Z. 

C 

COMPLEX  ZL(700)Z(S0000) 

NxMT+MP 
M0xMT3-2 
DO  100  1x1, MP 
L0xMrN+(l-1)*N*MT 
IF(I.EQ.MP)  GO  TO  80 
Klx2 

IF(I.E0.1)KI=1 

IF(I.EQ.MT)  Kl=3 

L2x(l-1)*N+I 

L1xL2-1 

L3xL2+1 

M=2+3*(l-2) 

IF(KI.EQ.1)  MxO 
MIxM-fl 
M2xM-f2 
M3xM>3 

go  to  (50.60,70)  .ki 
50Z(L2)xZ(L2)+ZL(M1) 


Z(L3)xZ(L3)«ZL(M2) 

GO  TO  80 

60Z(L1)xZ(L1)'i'ZL(M1) 

Z(L2)xZ(L2)'»ZL(M2) 

Z(L3)xZ(L3)*ZL(M3) 

GO  TO  80 

70Z(L1)>Z(L1)4'ZL(M1) 

Z(L2)xZ(L2)+ZL(M2) 

80  Z(L0+l)xZ(L0+l)+ZL(M0+l) 

100  CONTINUE 
RETURN 
END 

subroutin*  inl*(p(u,v,«.y.nn) 

c  program  to  inlatpolata  Inaarty  batwaan  arrays  of  x  arxi  y 
c  valuas  (with  difnansions  rm).  u  is  spacifiad  and  v  is 
ratumad 

c  unilonn  aamplii^  in  x  (tba  ordinata)  is  assutnad. 
dimansion  x(nn).y(nn) 
dalxxx(2)-x(1) 
n1xinl(u/dalx)'r'1 
nZxnl+l 
dalyxy(n2)-y(n1) 
sgnxsign(1..daly) 
dalyxabs(daly) 
a)phaxu-x(n1 ) 

2ataxdaly*alptwdalx 
vxy(n1  )«sgn‘zata 
ratum 
arrd 

SUBROUTINE 

PLANEW(NWIRES.NW1  .NW2,NP1  ,NP2.XH.YHZH, 
•THR.PHR.R) 

C 

C  PLANE  WAVE  EXCITATION  VECTOR  ELEMENTS  FOR 
WIRE  AND 

C  ATTACHMENT  SEGMENT.  INCIDENCE  DIRECTION  IS 
(THR,PHR). 

C  •  ONLY  ONE  INCIDENCE  ANGLE  PER  CALL 

C  -  WIRE  DOES  NOT  HAVE  TO  LIE  IN  THE  ZxO  PLANE 

C 

COMPLEX 

U0.AA,BB.C,R(400),CEXP.EXP.FI1,FI2.SI.DI.CMPLX 

DIMENSION 

ZH(400),NS(4).NW1(4).NW2(4),XH(400),YH(400) 
MP2=NP2-1 
MT2=NP2-2 
MT22»2'MT2 
D0  5L.1.NWIRES 
5  NS(L)xNW2(L)-NW1(L)+1 
U0x(0.,0.) 

CCxCOSfTHR) 

SSxSINfTHR) 

CPxCOS(PHR) 

SPxSIN(PHR) 

UPxSS'CP 

VPxSS’SP 

WPxCC 

DO  12  IPx1,MP2 
lIxIP+NPI 

1x114^1 

ZS..5'(ZH(I)+ZH(II)) 

XS>.S*(XH(I)+XH(II)) 

YS».S*(YH(I)*YH(II)) 

DXxXH(l)-XH(ll) 

DYxYH(l)-YH(ll) 
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01-SQRT(DX"2+DY**2) 

SU-DY/D1 

CU-DX/D1 

C  FOR  WIRES  IN  PERPENDICULAR  TO  Z  SIN(V).1  AND 
COS(V)»0 
SV-1.0 
CV.0.0 

WIRE  SEGMENT  CALCULATIONS 

A»UP’CU+VP’SU 
B»UP'XS+VP*YS+WP*ZS 
C»CMPLX(0  A) 

EXP»CEXP(CMPLX(0.,3)) 

AA=CC*(CU‘CP+SU*SP) 

BB»SU*CPSP‘CU 

PSI«D1*A/2. 

IF(PSI.NE.O.)  GO  TO  60 
SINC»1. 

GO  TO  61 

60  SINC>SIN(PSI)/PSI 

61  COSP«COS(PSI) 

FI1=SINC*DrEXP/2. 

FI2=(0.,0.) 

IF(ABS(A).LT.1.E-4)  GO  TO  62 
CSP-COSP-SINC 
IF(ABS(CSP).LT.1  .E-4)  GO  TO  62 
FI2=EXP/C'CSP 

62  CONTINUE 
SI=FIUFI2 
DI=FI1-FI2 

C 

C  R-WIRE-THETA 
C 

IF(IP  EO.MP2)  GO  TO  10 
B{IP)=AA'SI 
R(IP+MT2)=BB‘SI 
10  CONTINUE 
C 

C  R-WIRE-PHI 
C 

14IF(IP.EQ.1)GOT0  12 
R(IP-1).R(IP-1)+AA‘DI 
R{IP-1+MT2)»R(IP-1+MT2)+BB*OI 
12  CONTINUE 
C 

C  WRITE  OVER  DUMMY  WIRE  SEGMENTS 
C 

c  IF(NWIRES.EQ.I)  GO  TO  210 
c  NSPTS=0 
c  DO209L»1.NWIRES-1 
c  NSPTS=NSPTS+NS(L) 
c  IN=NSPTS-2 
c  R(IN+1).U0 
c  R(IN+2)sU0 
c  IN«IN+MT2 
c  R(IN+1)sU0 
c  R(IN+2)>:U0 
c  209  CONTINUE 
210  RETURN 
END 

SUBROUTINE 

PLANES(M1  .NPO.NPI  ,NP2,NT.RH2H,XTAT,THR,R) 

C 


C  PLANE  WAVE  EXCITATION  VECTOR  FOR  BOR 
ELEMENTS 

C  NO  CHANGE  FROM  HARRINGTON'S  ORIGINAL 
PROGRAM  OTHER  THAN 
C 

C  •'  ONLY  ONE  MODE  PER  CALL  " 

C 

COMPLEX 

R(1000),U,U1,UA.UB.FA(50).FB(50).F2A,F26.F1A.F1B 
COMPLEX  U2.U3.U4.U5.CMPU 
DIMENSION  RH(400)2H(400).XT(4)AT(4) 
DIMENSION  R2(4)22(4).BJ(5000) 

U>(0.1.) 

U1>3.141593‘U*‘M1 

NF<1 

M2>M1 

NP«NP0+NPUNP2 

NRF.NP0+NP1 

MP»NP-1 

MT-MP-1 

NmMT+MP 

N2»2-N 

CC=COS(THR) 

SS=SIN(THR) 

M3>MU1 

M4>M2+3 

IF(MI.EQ.O)  M3=2 

M5<MU2 

M6>M2+2 

D012IP«1.MP 

K2«IP 

lalP-fl 

DR..S-<RH(I).RH(IP)) 

DZ-.5*(ZH(I)-ZH(IP)) 

D1=SQRT(DR'DR+DZ*DZ) 

R1..25*(RH(I)+RH(IP)) 

IF(ABS(R1).LT.1.E-5)R1»1. 

Z1..S'(ZH(I)*ZH(IP)) 

DR..5*DR 
D2=DR/R1 
DO  13L»1.NT 
R2(L)=R1+DR*XT(L) 

Z2(L)»Z1+DZ*XT(L) 

13  CONTINUE 
D3=DR’CC 
D4.-DZ*SS 
D5«DrCC 
DO  23  M>M3.M4 
FA(M)»0. 

FB(M).0. 

23  CONTINUE 
D015L»1,NT 
X«SS‘R2(L) 

IF(X.GT..5E-7)  GO  TO  19 
DO  20  M>M3.M4 
BJ(M)>0. 

20  CONTINUE 
BJ(2)-1. 

S»1. 

GOTO  18 
19  Mi^,8'X+14.-2./X 
IF(X.LT..5)  M»11.8+ALOG10(X) 

IF(M.LT.M4)  M>M4 
BJ(M)s0. 

JM>M-1 
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DO  16  J=4  .M 
J2>JM 

BJ(JM)=J1/X‘BJ(J2)-BJ(JM+2) 

16  CONTINUE 
SaO 

IF(M.LE.4)  GO  TO  24 
DO  17  J=4,M,2 
S>S-fBJ(J) 

17  CONTINUE 

24  S>BJ(2)'»2.'S 

18  ARG=Z2(L)”CC 

UA=AT(L)/S‘CMPLX(COS(ARG),SIN(ARG)) 
UB»XT(L)'UA 
DO  25  MsM3.M4 
FA(M)=BJ(M)'UA+FA(M) 

FB(M)sBJ(M)‘UB>FB(M) 

25  CONTINUE 
15  CONTINUE 

IF(MI  .NE.O)  GO  TO  26 
FA(1)=-  FA(3) 

FB(1)=-  FB(3) 

26  UA=U1 

DO  27  M=M5.M6 

M7=M-1 

M8=M+1 

F2A=UA'(FA(M8)+FA(M7)) 

F2B=UA-(FB(M8)+FBlM7)) 

UB=U-UA 

F1A=UB-(FA(M8)-FA(M7)) 

F1B=UB'(FB(M8)-FB(M7)) 

U4=D4‘UA 

U2.D3'F1A+U4'FA(M) 

U3=D3’F1B*U4‘FB(M) 

U4.DRT2A 

U5=OR’F2B 

K1=K2-1 

K4=K1+N 

K5«K2+N 

R(K2+MT)=-D5'(F2A+D2'F2B) 

R(K5+MT)=Dr(F1A+D2’F1B) 

IF(IP.EQ.1)G0T0  21 
R(K1)=R(K1)+U2-U3 
R(K4)=R(K4)+U4-U5 
IF(IP.EQ.MP)  GO  TO  22 

21  R(K2)=U2+U3 
R(K5)=U4+U5 

22  K2=K2+N2 
UA=UB 

27  CONTINUE 
12  CONTINUE 

c  R(NP0-1)=(0..0.) 
c  R(NP0)=(0..0.) 
c  R(NP0+MT)=(0..0.) 
c  R(NRF-1)»(0..0.) 
c  R(NRF)='0..0.) 
e  R(NRF+MT)=(0..0.) 

99  RETURN 
END 

SUBROUTINE 

ZMATSS(M1  .M2.NP0.NP1  .NP2.NPHI.NT,RH^H.X.A.XT.A 
TZ) 

C 


C  ••••”*•  THREE  DETACHED  SURFACES  . 

C  SURFACE-SURFACE  IMPEDANCE  ELEMENTS 
REMAINS  UNCHANGED  FROM  HARRINGTON 
C  EXCEPT  multiple  SURFACES  ARE  PERMITTED 
THE  FIRST  SURFACE  HAS  NPO 
C  POINTS.  THE  SECOND  NP1  POINTS  AND  THE  THIRD 
NP2 

C  NPO  PARABOLOID  PTS:  NP1  HYPERBOLOID  PTS, 

NP2  CAVITY  PTS. 

C 

COMPLEX 

Z(50000).U1.U2.U3.U4.U5.U6.U7.U8.U9.UA.UB.G4A(4|  G5 
A(4) 

COMPLEX 

G6A(4).G4B(4).G5B(4).G6B|4).H4A.H5A.H6A.H4B.HSB 
COMPLEX  CMPLX.H6B.UC.UO.GA|400|.GB(400| 
DIMENSION 

RH(400)ZH(400)  .X(400).A(400) .  AT(4)  .RS(400)  ZS|400| 
DIMENSION 

D(400).DR(400).OZ(400).DM(400).C2(400).C3(400).R2(4) 

DIMENSION 

C4(400)  .05(400)  .06(400)  Z7(4)  .R7(4)  Z8(4)  .R8(4) 
DIMENSION  XT(4)Z2(4) 

CT.2. 

CP=  1 

NP=NP0+NP1+NP2 
NRF»NP0-fNP1 
DO  10  l=2.NP 
I2sl-1 

RS(I2)«.5"(RH(I)»RH(I2)) 

ZS(I2)..5*(ZH(I)+ZH(I2)) 

01  =  .5‘(RH(I)-RH(I2)) 

D2..5’(ZH(I)-ZH(I2)) 

D(l2)»SORT(D1  ‘D1  +D2‘D2) 

DR(I2).D1 

0Z(I2)>02 

0M(I2).D(I2)'RS(I2) 

10  CONTINUE 
M3<M2-MU1 
M4sM1-1 
PI2. 1.570796 
DO  11  K»1,NPHI 
PH«PI2’(X(K)+1.) 

C2(K)=PH'PH 

SN=SIN(.5*PH) 

C3(K).4.’SN’SN 

A1=PI2-A(K) 

D4..5'A1'C3(K) 

D5.ArCOS(PH) 

D6-A1'SIN(PH) 

M5>K 

D0  29M-1.M3 

PHM«(M4+M)’PH 

A2-COS(PHM) 

C4(M5)>I34*A2 

C5(M5)-D5’A2 

C6(M5)>D6*SIN(PHM) 

M5sM5-fNPHI 
29  CONTINUE 

11  CONTINUE 
MP=NP-1 
MT.MP-1 
N=MT+MP 
N2N.MT’N 
N2=N’N 
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U1«(0.,.u/ 

U2-(0.,2.) 

JN--1-N 
DO  15  JQ>1.MP 
KO-2 

IF(JQ.EO.I)  KQ«1 
IF(JQ.EQ.MP)  KQ»3 
R1>RS(ja) 

Z1>ZS(JQ) 

01«0(JQ) 

D2>OR(JQ) 

D3>0Z(JQ) 

D4>02/R1 

DS>DM(JQ) 

SV-D2/D1 

CV-D3/D1 

T6«croi 

T62-T6+D1 

Tsz-TezTez 

R6-CP-R1 

R62«R6'R6 

DO  12L»1,NT 

R2(L)-R1+D2-XT(L) 

Z2(L)=Z1+D3-XT(L) 

12  CONTINUE 
U3=D2-U1 
U4-D3'U1 
DO  16  IP=1,MP 
R3>RS(IP) 

Z3^S(IP) 

R4.R1-R3 

Z4=Z1-Z3 

FM«R4'SV+Z4‘CV 

PHM>ABS(FM) 

PH-ABS(R4*CV-Z4'SV) 

D6.PH 

IF{PHM,LE.D1)GOT0  26 
D6=PHM-01 

D6-SQRT(D6’D6+PH-PH) 

26  IF(IP.EQ.JQ)  GO  TO  27 
KP=1 

IF(T6.GT.D6)  KP=2 
IF(R6.GT.D6)  KP=3 
GO  TO  28 

27  KP=4 

28GOTO{41,42.41.42).KP 

42DO40L=1,NT 

D7=R2(L)-R3 

08>Z2(L)-Z3 

Z7(L)=D7'D7+D8’D8 

R7{L)=R3*R2(L) 

Z8{L)=.25‘Z7(L) 

R8(L)«.25*R7(L) 

40  CONTINUE 
Z4=R4‘R4+Z4*Z4 
R4«R3’R1 
R5-.5‘R3*SV 
D0  33K»1,NPHI 
A1>C3(K) 

RR>Z44^R4‘A1 

UA-0. 

UB>0. 

IF(RR.LT.T62)  GO  TO  34 
D0  3SL»1,NT 
R»SQRT(Z7(L)+R7(L)*A1) 


SN>-SIN(R) 

CS>COS(R) 

UC«AT<L)/R-CMPU(CS,SN) 

UA«UA«UC 

UB>)a(L)'UC+UB 

35  CONTINUE 
GO  TO  36 

34D0  37L-1,NT 
R«SQRT(Z8(L)+R8(L)*A1 ) 
SN>-SIN(R) 

CS«COS(R) 

UC«AT(L)/R*SN‘CMPU(-SN.CS) 

UAsUA+UC 

UB«)n’(L)‘UOUB 

37  CONTINUE 
A2>FM+R5'A1 
D9>RRA2‘A2 
R>ABS(A2) 

D7-R-D1 

D8sR+D1 

06>SQRT(D8'D8+D9) 

R. SQRT{D7-07+D9) 

|F(D7.GE.0.)  GO  TO  38 

AI. AUOG((D8+D6)*(-D7+R)/08)'D1 
GO  TO  39 

38  A1>ALOG((D84'D6)/(D7'fR))/D1 

39  UA^AUUA 

UB>A2'(4./(06«R)-A1  )/DUUB 

36  GA(K).UA 
GB<K)-UB 

33  CONTINUE 
K1«0 

D0  45M>1.M3 
H4AxO. 

H5A>0. 

H6A>0. 

H4BsO. 

HSB>0. 

H6B>0. 

D0  46K.1.NPHI 

K1=KU1 

D6xC4(K1) 

D7.C5(K1) 

De=C6(K1) 

UA=GA(K) 

UB«GB(K) 

H4AsD6’UA+H4A 
H5A.07*UA+H5A 
H6A«D8*UA+H6A 
H4B«D6'UB-»H4B 
H5B.D7’UB+H5B 
H6B«D8'U&fH6B 
46  CONTINUE 
G4A(M)=H4A 
G5A(M)sH5A 
G6A(M)>H6A 
G4B(M)>H4B 
G5B(M)irfH6B 
G6B(M)«H6B 
45  CONTINUE 
|F(KP.NE.4)  GO  TO  47 
A2»D1/(PI2'R1) 

D6-2./D1 

D8>0 

D0  63K-1,NPHI 
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A1«R4'C2(K) 

R>R4'C3(K) 

IF(R.LTT62)  GO  TO  64 
D7»0. 

D0  65L-1.NT 

D7.D7+AT(L)/SQRT(27(L)+A1 ) 

65  CONTINUE 
GO  TO  66 

64  A1bA2/(X(K)+1.) 

D7=D6*ALOG(A1  ♦SQRT{  1  .+A1  ‘A1 )) 

66  D8»D8+A(K)’D7 
63  CONTINUE 

A1«.S’A2 

A2»1./A1 

D8=-PI2’D8*2./R  1  '(BLOG{A2)+A2*BLOG(A1 )) 

00  67  M«1,M3 
GSA(M)«D8-fGSA(M) 

67  CONTINUE 
GO  TO  47 

41  D0  2SMs1,M3 
G4A(M)=0. 

G5A(M)=0. 

G6A(M)sO. 

G4B(M)sO. 

G5B(M)=0. 

G6B(M)=0. 

25  CONTINUE 
DO  13  L=1.NT 
A1=R2(L) 

R4=A1-R3 

24=Z2(L)-Z3 

24»R4'R4+24*24 

R4=R3*A1 

DO  17K=1,NPHI 

RsSQRT(24+R4"C3(K)) 

SN=-SIN(R) 

CS«COS(R) 

GA(K)=CMPLX(CS.SN)'R 
17  CONTINUE 
D6=0. 

IF(R62.LEi4)  GO  TO  51 
DO  62  K=1,NPHI 

D6=D6+A(K)/SQRT(24+R4*C2(K)) 

62  CONTINUE 
24=3. 1 41 593/SQRT(24/R4) 

D6=-PI2*D6+ALOG(24+SQRT(  1  .+24‘24))/SQRT(R4) 
51  A1=AT(L) 

A2=XT(L)‘A1 

K1=0 

DO  30  M=1,M3 
U5=0. 

U6=0. 

U7=0, 

D0  32K=1,NPHI 
UA=GA(K) 

K1.K1+1 

U5»C4(K1)*UA+U5 
U6>C5(K1)*UA+U6 
U7»C6(K1)*UA+U7 
32  CONTINUE 
U6=06+U6 

G4A(M)=A1  'U5+G4A(M) 

G5A(M)=A1*U&fG5A(M) 

G6A(M)-A1  ■U7+G6A(M) 

G4B(M)>A2*U&fG4B(M) 


G5B(M)>A2’Ue+G5B(M) 

G6B(M)=A2'U7+G6B(M) 

30  CONTINUE 

13  CONTINUE 

47A1«DR(IP) 

UA=ArU3 

UB>D2(IP)*U4 

A2>D(IP) 

D6.A2"02 

D7«DrA1 

D8>DrA2 

JM>JN 

DO  31  M=1,M3 

FM>M4+M 

A1.FM'DM(IP) 

H5A>G5A(M) 

H5B>GSB(M) 

H4A=G4A(M)+H5A 

H4B=G4B(M)+H5B 

H6A=G6A(M) 

H6B=G6B(M) 

U7=UA*HSA+UB‘H4A 

U8=UA*H5B«UB‘H4B 

U5=U7-U8 

U6<U7+U8 

U7.-UrH4A 

U8=D6’H6A 

U9=D6*H6B-ArH4A 

UC=D7‘(H6A4'D4*H6B) 

UD=FM*D5‘H4A 

KUIP+JM 

K2=K1+1 

K3=K1+N 

K4cK2+N 

K5=K2+MT 

K6-K4+MT 

K7=K3+N2N 

K8=K44.N2N 

GOTO(18,20.19).KQ 

18  2(K6)=U8+U9 
IF(IP.EQ.1)G0T0  21 
2(K3)=2(K3)+U6-U7 
2(K7)=2(K7)+UC-UD 
IF(IP,EQ.MP)  GO  TO  22 

21  2(K4)=U6+U7 
2(K8)=UC+UD 
GO  TO  22 

19  2(K5)=2(K5)+U8-U9 
IF(IP.EQ.1)G0T0  23 
2(K1)=Z(K1)+U5+U7 
2(K7)=2(K7)+UC-UD 
IF(IP.EQ.MP)  GO  TO  22 

23  2(K2)=2(K2)+U5-U7 
2(K8)=UC+UD 

GO  TO  22 

20  2(K5)>2(K5)i'U8-U9 
2(K6)=U84^U9 
IF(IP.EQ.1)G0T0  24 
2(K1).2(K1)+U5+U7 
2(K3)=2(K3)+U6-U7 
2<K7)«2(K7)+UC-UD 
IF(IP.EQ.MP)  GO  TO  22 

24  2(K2)»2(K2)+U5-U7 
2(K4)=U6+U7 
2(K8)>UC■^UD 
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22  Z(K8+MT)«U2*(D8'<H5A+D4*H5B)-A1*UD) 
JMsJM+N2 
31  CONTINUE 
16  CONTINUE 
JNaJN-fN 
IS  CONTINUE 
C 

C  THREE  MULTIPLE  BODIES  USING  THE  SIMPLIFIED 
APPROACH 

C  NULL  OUT  ROWS  AND  COLS  FOR  THE  FIRST 
DUMMY  SEGMENT 
C 

DO  100LSS-1.N 
LSsLSS-1 

Z(LS*N+NP0)=(0..0.) 

Z(LS*N+NPO-1)=(0.,0.) 

100  Z(LS*N+NPO+MT)=(0.,0  ) 

DO  101  LS=1.N 
Z((NP0-2)’N+LS)=(0..0.) 

Z((NPO-1)’N+LS)=(0..0  ) 

101  Z((NP0-1+MT)*N+LS)=(0..0.) 

Z((NPO-2)‘N+NPO- 1  )=(  1  .0.) 
Z((NPO-1)'N+NP0)=(1..O  ) 

Z((MT+NP0- 1  )*N+NP0+MT)=(  1 .0.) 

C 

C  NULL  OUT  ROWS  AND  COLS  FOR  THE  SECOND 
DUMMY  SEGMENT 
C 

DO  200  LSS=1,N 
LS=LSS-1 

Z(LS*N+NRF)=(0  .0  ) 

Z(LS-N+NRF-1)»(0..0) 

200  2(LS*N+NRF+MT)=(0..0.) 

DO  201  LSsl.N 
2((NRF-2)'N+LS)=(0..0.) 

Z((NRF-1)’N+LS)=(0..0.) 

201  Z((NRF-1+MT)-N+LS)*(0..0.) 
Z((NRF-2)*N+NRF-1)=(1..0,) 
2((NRF-1)*N+NRF)=(1..0,) 

2((MT+NRF- 1  )'N+NRF+MT)=(  1 .0  ) 

999  RETURN 
END 

SUBROUTINE  ZASMB0(NP1  .NP2.MODES^ZWW) 
C 

C  ASSEMBLES  MODE  INDEPENDENT  BLOCKS  OF  Z 
C 

COMPLEX  Z(500000)ZWW(10000) 

MT1=NP1-2 
MP1=NP1-1 
N=MT1+MP1 
MT2=NP2-2 
NBLOCK=2*MODES+1 
NS=NBLOCK-N 
NR=NS+MT2 
NS2=NS*NR 
NWW=NS2+NS 
DO  10  IC«1.MT2 
IC1=IC-1 
DO  10  IR=1,MT2 
K=NWW+ICrNR+lR 
IW«ICrMT2+IR 
10  Z(K)«2WW(IW) 

RETURN 

END 


SUBROUTINE 

ZASMBN(NP1  ,NP2.MODES.NMZi:sS^SW) 

C 

C  ASSEMBLES  MODE  DEPENDENT  BLOCKS  IN  Z 
MODES>TL  •  OF  MODES; 

C  NM  IS  THE  MODE  •  OF  THE  BLOCKS  ZSS2SW2SJ 
(IF  NM<0  THEN 

C  ZSS  IS  OMITTED  •  ONLY  FILLED  FOR  POSITIVE  NM) 
C 

COMPLEX  ZSS(SOOOO);!(50OOOO|^SW(50OOO) 

MT1.NP1-2 

MP1=NP1-1 

MT2.NP2-2 

N>MP1+MT1 

NBLOCK=2'MODES+1 

NS>NBLOCK'N 

NRsNS'fMTZ 

MT-MT1+MT2 

MTN=MT*N 

IF(NM.LT.O)  GO  TO  225 
C 

C  ZSS  BLOCKS  FOR  MODE  NM  (+  AND  -  FILLED  IN  1 
LOOP) 

C 

NTOP1=MODES-NM 
NTOP2=NBLOCK(NTOP  1  ♦  1 ) 

NSS2=NTOP1  ‘(N+NR'N) 

NSS1=NTOP2'(N+NR*N) 

DO  200  IC.1.MP1 
IC1=IC-1 

DO  200  IR=1,MP1 
IF(IC.GT.MT1)GOT0  100 
IF(IR.GT.MT1)GOTO50 
IO=ICrN+IR 

rrr-''^  'nr+ir 

Z(r.*NSSi)*zss((0) 

IF(NM.NE.0)  Z(rrT+NSS2)=ZSS(IO) 

50  I0=MT1+IC1'N+IR 
IPT=ICrNR+MTUIR 
Z(IPT+NSS1)=  ZSS(IO) 

IF(NM.NE.0)  Z(IPT+NSS2)=-ZSS(I0) 

100  IFdR  GT.MT1)  GO  TO  150 
IO=N‘MTUICrN+IR 
ITPxMTI  ’NR+ICI  "NR+IR 
Z(ITP+NSS1)=  ZSS(IO) 

IF(NM.NE.O)  Z(rTP+NSS2)=-ZSS(l0) 

150  IO=N'MT1+MT1+ICrN+IR 
IPP=ICVNR+MT1+MTrNR+IR 
Z(IPP+NSS1)=ZSS(I0) 

IF(NM.NE.O)  Z(IPP+NSS2)»ZSS(I0) 

200  CONTINUE 
225  CONTINUE 
C 

C  ZSW  AND  ZWS  BLOCKS 
C 

NTOP1«MODES+NM 

NTOP2«NBLOCK-(NTOP1+1) 

NSW.NTOP  1  ’N+NSHR 

NWS«NTOP2'NR’N+NS 

DO250IC«1,MT2 

ICUIC-1 

DO250IR-1.N 

ISW.IR+ICVNR+NSW 

IWS«(IR-1)’NR+IC+NWS 

lO^ICl’N+IR 
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Z(ISW)-ZSW(IO) 

Z(IWS)>ZSW(IO) 

2S0  CONTINUE 
RETURN 
END 

SUBROUTINE 

ZMATSW(NWIRES,NW1  ,NW2.NP1 1  ,NP12,NP13, 
•XH,YH,RHZH,NT.XT^T,NGJ<G.AG,MODE.AZ) 

C 

c  »»»»>»  r*vB*d  version  ««««« 
c 

C  CALCULATION  OF  THE  SURFACE-WIRE  MATRIX 
ELEMENTS.  ATTACHMENT 

C  POINT  NOT  AT  A  SURFACE  EDGE,  BUT  CAN  BE  ON  Z 
AXIS  (IC«1).  SURFACE 

C  START  INDEX  IS  1 .  TWO  BORS  OF  LENGTHS  NP1 1 , 
NP12AND  NP13.  WIRE 

C  START  AND  STOP  POINTS  ARE  NW1  AND  NW2. 

C 

COMPLEX 

CEXP.EXPZ(50000),G5,G6.G7,CONH.F5.F6.F7 

COMPLEX 

U.UO.PSI,TO.T1,T2,T3,ST1.ST2.ST3.CON,SP1,SP2.SP3 
COMPLEX  U5,U6,U7.U8,U9.CMPLX,T4,TS,T6,T7 
DIMENSION 

RS1(400),RH(400)2H(400),XG(400),AG(400).XT(4),AT(4) 

DIMENSION 

ZS1(4CX)),D1(400),D(400).S(400),C(400),RS(400) 

DIMENSION 

ZS(400),S1(400),C1(400),NW1(4).NW2(4),NS(4) 

DIMENSION 

XH(400),YH(400)  .XS1  (400)  ,YS  1  (400)  .CU(400)  ,SU(400) 
Pl=3,14159 
PI2=2.”PI 

NP1aNP11+NP12+NPt3 
NRF=NP11+NP12 
DO  1  L»1.NWIRES 
1  NS(L)=NW2(L)-NW1(L)+1 
U0=(0..0.) 

U=(0..1.) 

C 

C  DEFINE  SURFACE  GEOMETRY  TERMS 
C 

MP=NP1-1 
MT=MP-1 
NROW=MP+MT 
DO  5  l=2.NP1 
12=1-1 

RS(I2)=.5*(RH(I)+RH(I2)) 

ZS(I2)=.S*(ZH(I)+ZH(I2)) 

DR=RH(I)-RH(I2) 

DZ=ZH(I)-ZH(I2) 

D(I2)=SQRT(DR"2+DZ"2) 

S{I2)=DR/D(I2) 

5  C(I2)=DZ/D(I2) 

C 

C  DEFINE  WIRE  GEOMETRY  TERMS 
C 

NW=NW2(NWIRES)-NW1  ( 1 ) 

NWP=NW+1 

MT2=NWP-2 

DO  10  N»2.NWP 

N0=N-1 

l=NP1'fN 

12=1-1 


RS1(N0)«.5‘(RH(I)4'RH(I2)) 

ZS1(N0)».5*(ZH{I)+ZH(I2)) 

XS1(N0)>.5‘(XH(I)+XH(I2)) 

YS1(N0)..5'(YH(I)+YH(I2)) 

DX>XH(I)-XH(I2) 

DY«YH(I)-YH(I2) 

D1(N0)=SQRT{DX"2+DY"2) 

C  FOR  WIRES  PARALLEL  TO  THE  XY  PLANE  SIN(V).1 
AND  COS(V)«0 
S1(N0)>1. 

C1(NO)>0. 

CU(N0)>DX/D1(N0) 

SU(N0)>DY/D1(N0) 

10  CONTINUE 
CON«U/PI2(2. 

C  INTEGRATING  FROM  0  TO  2‘PI 
01«PI 
02»PI 

JN-1-NROW 

C 

C  SEGMENT  LOOPS:  q=WIRE.  p=SURFACE 
C 

DO  100  JQ=1.NW 
K(3=2 

IF(JQ.EQ.1)  KQ=1 
IF(JQ.EQ.NW)  KQ=3 
H1=D1(JQ)/2. 

ZW=ZS1(JQ) 

D0  22IP=1.MP 

ST1.U0 

ST2=U0 

ST3=U0 

SP1=U0 

SP2=U0 

SP3=U0 

AA=S(IP)-CU(JQ) 

BB=S(IP)*SU(JQ) 

RD=RS(IP)'D1(JQ) 

0D«D(IP)*D1(JQ) 

CONH=CON’D(IP) 

C 

C  FIRST  TERM  (H-INTEGRATION) 

C 

DO  30  1.1, NT 
H.HrXT(l) 

HHD.2.'H/D1(JQ) 

XW.XS1(JQ)+H*CU(JQ) 

YW=YS1(JQ)+H*SU(JQ) 

C 

C  PHI  INTEGRATIONS  FOR  G-HAT  FUNCTIONS 
C 

G7=U0 

G5=U0 

G6.U0 

D0  32K.1,NG 

PH.OrXG(K)+02 

PHM.PH*FLOAT(MODE) 

XS«RS(IP)"COS(PH) 

YS.RS(IP)'SIN(PH) 

EXP.CEXP(CMPU(0„-PHM)) 

RP.SQRT((XS-XW)"2+(YS-YW)**2+(ZS(IP)-ZW)“2) 

PSI.CEXP(CMPLX(0.,-RP))/RP 

F5«EXP*COS(PH)*PSI 

F6=EXP*SIN(PH)*PSI 

F7-EXP*PSI 
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o  o  o  o  o 


G7.G7+F7'AG(K) 

G5>GS+FS'AG(K) 

32  G6«GfrfF6'AG(K) 

G7«OrG7 

G&iOI'GS 

G6-OrG6 

T0«GS*AA/4. 

TI-TO'HHD 

T2«G6’BB/4. 

T3«T2*HHD 

T4»G5-SU(JQ)/2. 

T5»T4‘HHD 

T6.-G6*CU(JQ)/2. 

T7.T6*HHD 

ST1  .ST1  ♦AT(I)*(T0+T2) 

ST2-ST2+AT(1)*(TUT3) 

ST3-ST3+AT(I)"G7/DD 

SP1.SPUAT(I)*(T4+T6) 

SP2»SP2+AT(I)*(T5+T7) 

SP3-SP3+AT(l)*U*G7*FLOAT(MODE)/RD 

30  CONTINUE 
C 

C  ZSW-t  TERMS  (K1  ,K2,K3.K4) 

U5=(ST1  -ST2)*CONH‘H  1 
U6=(ST1 +ST2)"CONH‘H  1 
U7=-ST3*CONH-H1 

ZSW-p  TERMS  (KS.K6) 

U8=SPrCONH'H1 

U9=(SP2+SP3)*CONH’H1 

KUIP+JN 

K2*KU1 

K3aK  UNROW 

K4*K2+NROW 

KSSK2+MT 

K6«K4+MT 

GOTO(18.20.19).KQ 

18Z(K6)=U8+U9 
IF{IP.EQ.1)G0T0  21 
Z(K3)=Z(K3)+U6-U7 
IF(IP.EQ.MP)  GO  TO  22 

21  Z(K4)«U6+U7 
GO  TO  22 

19Z(K5)=Z(K5)+U8-U9 
IF(IP.EQ.1)G0T0  23 
Z(K1)=Z(K1)+US+U7 
IF(IP.EQ.MP)  GO  TO  22 

23  Z(K2)=Z(K2)+U5-U7 
GO  TO  22 

20Z(K5)<2(K5)+U8-U9 
Z(K6)sU8-fU9 
IF(1P.EQ.1)G0T0  24 
Z(K1)=Z(K1)+US+U7 
Z(K3)=Z(K3)+U6-U7 
IF(IP.EQ.MP)  GO  TO  22 

24  Z(K2)=Z(K2)+U5-U7 
Z(K4)«U6+U7 

22  CONTINUE 
JN>JN-)-NROW 

100  CONTINUE 

FOR  MULTIPLE  BORS  WRrTE  OVER  DUMMY 
SEGMENT  ROWS 


C  VmrTE  OVER  ROWS  AND  COLS  FOR  FIRST  DUMMY 

SEG 

C 

DO  101  LSS«1,MT2 
LS-LSS-1 

ZILS-NROW+NP1 1)«(0.,0.) 

Z(LS"NROW+NP1 1-1  )«(0.,0.) 

101  Z(LS'NROW+NP11-fMT)>(0..0.) 

C 

C  VI/RITE  OVER  ROWS  AND  COLS  FOR  SECOND 
DUMMY  SEG. 

C 

DO  201  LSS>1.MT2 
LS«LSS-1 

Z(LS‘NROW+NRF).(0..0.) 

Z(LS*NR0W4-NRF-  1  )>(0..0.) 

201  Z(LS‘NROW+NRF+MT)>(0..0.) 

C 

C  WRITE  OVER  DUMMY  SEGMENTS  IN  THE  CASE  OF 
MULTIPE  WIRES 
C 

IF(NWIRES.EQ.  1 )  GO  TO  292 
NSPTS=0 

DO  291  L=1.NWIRES-1 
NSPTS=NSPTS+NS(L) 

DO  291  l>1.NROW 
C 

C  SET  COLS  TO  ZERO  FOR  STRIPS 
C 

IN=(NSPTS-2)*NROW 

Z(IN+l)>U0 

Z(IN-t.|+NROW)>U0 

291  CONTINUE 

292  RETURN 
END 

SUBROUTINE 

ZMATWW(NWIRES,NW1,NW2,XH.YH.RH2H,NT.XT.AT,A, 

2) 

C 

C  •••  MODS  FOR  DIPOLE  --  RECTANGULAR 
COORDINATES  ARE  SENT  OVER  **• 

C  COMPUTE  MATRIX  ELEMENTS  FOR  MULTIPLE 
WIRES  USING  GLISSON  METHOD 
C  NWIRES-  #  OF  WIRES;  NW1(I).NW2(I)  ARE  START 
AND  STOP  POINTS. 

C  &&&&&  WIRES  NEED  NOT  BE  IN  THE  Z°0  PLANE 
S&&S& 

C 

COMPLEX  CEXPZ(10000),CON,CMPLX 
COMPLEX  U.UO.PSI.SUMI  ,SUM2,SUM3.US.U6.U7 
DIMENSION 

RH(400)ZH(400).XT(4),AT(4),XH(400),YH(400) 
DIMENSION  D1(400),S1(400).C1(400) 

DIMENSION 

NW1(4).NW2(4).NS(4).XS1(400).YS1(400)2S1(400) 
DIMENSION  CU(400).SU(400) 

Pte3.141592 

PI2=2.'PI 

PD2=PI/2. 

U0>(0..0.) 

U«<0.,1.) 

C 

C  DEFINE  GEOMETRY  TERMS  FOR  THE  WIRE. 
XH.YH,RH2H  ARE  ALL  KNOWNS. 

C 


65 


0O5L«1.NWIRES 
5  NS(L)«NW2(L)-NW1(L(+1 
NSlsNW2(NWIRES)-NW1  ( 1 ) 

NPS-NS1+1 
DO  10  N.2.NPS 
N0>N-1 
l»NW1(1)+N-1 
12-1-1 

XS1(N0)-.5*(XH(I)+XH(I2)) 

YS1(N0)».S'(YH(I)+YH(I2)) 

2S1(N0)>.5'(ZH(I)+ZH(I2)) 

DX»XH(I)-XH(I2) 

DY«YH(I)-YH(I2) 

D1  (N0)=SQRT(DX**2+DY”2) 

CU(N0)»DX/D1(N0) 

SU(NO)>OY/D1(NO) 

S1(N0)=1.0 
C1(N0)=0.0 
10  CONTINUE 
NROW-NS1-1 
JN—1-NROW 
DO500  JQ=1,NS1 
KQ=2 

IF(JQ,EQ.1)  KQ=1 

IF(JQ.EQ.NS1)KQ=3 

Q1=D1(JQ),2. 

DO  22  IP-1, NS1 
LQ-0 

IF(IP.EQ.JQ)  LQ=1 
CON=D1(IP)/4 /pru 
SUMl-UO 
SUU2-U0 
SUM3-U0 

AA«CU(IP)’CU(JQ)+SU(IP)'SU(JQ) 

C 

C  H  INTEGRATION  -  SUBTRACT  OUT  SINGULARITY  IF 
JQ-IP 

C  AND  SAME  SEGMENT  (N11-N21);  DESIGNATED  LQ-1 
C 

DO  100  1*1. NT 
H=QrXT(l) 

IF(LQ.NE  1)GOTO40 
RP=SQRT(H”2+A*A) 

GO  TO  50 

40  XP=XS1(JQ)+H*CU(JQ)'S1(JQ) 

YP*YS1  (JQ)+H*SU(  JQ)‘S1  (JO) 

RP-SQRT((XS1  (IP)-XP)"2+(YS1  (IP)-YP)"2+(ZS1  (JQ)-ZS 
1(IP))"2) 

50  PSI=CEXP(CMPLX(0..-RP)) 
e  IF(LQ.NE.1)GOTO60 
c  PSI=PSI-(1..0.) 

60  PSI-PSI/RP 
SUM1=SUM1+PSrAA*AT(l)/4. 
SUM2=SUM2+PSI*AA*AT(  l)/2./Dl  ( JQ)*H 
SUM3-SUM3*PSrAT(l)/D1  (IP).'DI  (JQ) 

100  CONTINUE 
X1-0. 

X2-0. 

X3-0. 

c  IF(LQ.NE.1)GOTO200 
c  DQ-D1(JQ)/2. 
c  SDQ-SQRT(DQ**2+A**2) 
c  XIO-ALOG(SOO+Oa)-ALOG(SDa-DQ| 
c  X1*XI0*AA/4, 


e  X3-XI0.O1(IP).t)1(JQ) 
c  200  CONTINUE 

SUM1-(SUM1'QUX1)'CON 

SUM2«(SUM2'QUX2)’CON 

SUM3-(SUM3'Q1«X3)‘CON 

US-SUM1-SUM2 

U6-SUMUSUM2 

U7-SUM3 

K1-IP+JN 

K2-K1*1 

K3-K  UNROW 

K4«K2*NROW 

GOTO  (18.20, 19).KQ 

18  CONTINUE 
1F(IP.EQ  1)GOT0  21 
Z(K3)«Z(K3)«U6-U7 
IF(IP.EQ.NS1)G0T0  22 

21  Z(K4).U6*U7 
GO  TO  22 

19  CONTINUE 
IF(IP.EQ.1)G0T0  23 
Z(K1)-Z(K1)+U5+U7 
IF(IP.EQ.NS1)GO  TO  22 

23Z(K2)=Z(K2)*U5-U7 
GO  TO  22 

20  CONTINUE 
IF(IP.EQ  1)GO  TO  24 
Z(K1)»Z(K1)+U5fU7 
Z(K3)-Z(K3)+U6-U7 
IF(IP.EQ.NSI)  GO  TO  22 

24  Z(K2)»Z(K2)+U5-U7 
Z(K4)=U6+U7 

22  CONTINUE 
JN-JN+NROW 

500  CONTINUE 
IF(NWIRES.EQ.I)  GO  TO  297 
NSPTS-0 

DO  291  L*1,NWIRES-1 
NSPTS«NSPTS+NS(L) 

DO  291  I-1.NROW 
C 

C  USE  THE  SIMPLIFIED  METHOD  OF  HANDLING 
MULTIPLE  WIRES 

C  SET  COLS  TO  ZERO  FOR  STRIPS 
C 

IN.(NSPTS-2)’NROW 

Z(IN+l)-U0 

Z(IN+I+NROW)-U0 

C 

C  SET  ROWS  TO  ZERO  FOR  STRIPS 
C 

IN=NSPTS-2+(  I- 1  )*NROW 
Z(IN4^1)*U0 
Z(IN+2)-U0 
291  CONTINUE 
C 

C  SET  DIAGONALS  TO  1  FOR  STRIPS 
C 

NSPTS-0 

D0  296L-1.NWIRES-1 
NSPTS-NSPTS+NS(L) 

10- NSPTS-2 

11- IO*NROW+NSPTS-2 
ZdUlHT.O.) 

Z<IUNROW+2)=(1.,0.) 
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296  CO^mNUE 

297  RETURN 
END 
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